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Final  Report  on  ONR  Contract  N00014-75-C-0537 
Title:  Gulf  Oceanography-Gulf  Stream  Study  (GUSS) 

Project  3200A-6 

Principal  Investigator:  A.  D.  Kirwan,  Jr.,  Texas  A&M  University 

I.  Introduction 

In  1975  the  Office  of  Naval  Research  initiated  a theoretical 
and  observational  study  of  the  dynamics  of  the  Gulf  Stream.  Two 
phases  were  planned  for  the  study.  The  first,  or  feasibility  phase, 
was  to  be  conducted  during  1975-1978  with  the  observational  phase 
commencing  in  1979.  The  first  phase  has  a technical  and  scientific 
component.  The  former  is  being  conducted  by  the  Charles  Stark 
Draper  Laboratory,  Inc.  under  the  direction  of  W.  A.  Vachon.  The 
scientific  component  is  the  responsibility  of  A.  D.  Kirwan,  Jr.  of 
TAMU. 

Because  of  the  successes  reported  by  Kirwan  £t_  al_. , (1976)  and 
Richardson  (1976)  in  tracking  the  Culf  Stream  and  Gulf  Stream  rings, 
it  was  felt  that  the  primary  observational  tool  of  the  study  would 
be  Lagrangian  measurements  of  trajectories,  velocities,  accelerations, 
differential  kinematic  properties  (DKP)  such  as  horizontal  divergence 
and  vorticity  and  Reynolds  stresses.  During  the  first  year  of  the 
feasibility  phase,  the  A&M  group  performed  an  analysis  of  the  errors 
accruing  in  calculations  of  the  above  quantities  from  random  position 
error  of  Lagrangian  drifters  and  the  rate  from  which  the  position  was 
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determined.  In  addition,  some  preliminary  field  work  was  undertaken. 
The  Draper  group  performed  a study  of  position  fixing  systems  to 
determine  which,  if  any,  could  meet  the  accuracy  goals  as  specified 
by  the  analyses,  and  yet  still  be  cost  effective.  This  is  a final 
report  of  the  first  years  activities  on  this  project. 

The  details  of  the  preliminary  field  studies,  error  analysis, 
and  position  fixing  study  are  given  in  appendices  A,  B,  and  C, 
respectively.  Part  of  the  research  described  in  appendices  A and 
B was  supported  by  the  Office  for  the  International  Decade  of  Ocean 
Exploration  of  the  National  Science  Foundation. 

The  next  section  summarizes  the  findings  of  each  of  these 
studies  and  presents  our  recommendations  for  future  developments. 

II . Summary  of  Conclusions  and  Recommendations 
A.  Summary  of  conclusions 

1)  The  pilot  experiment  in  the  Kuroshio  described  in  Appendix  A 
demonstrated  that  the  drifters  stay  in  the  region  of  strong 
flowing  currents  long  enough  to  determine  all  of  the  required 
kinematic  parameters.  The  trajectories  and  velocities  ob- 
tained from  the  drifters  are  in  excellent  agreement  with  the 
hydrography. 

2)  Appendices  B and  C show  that  reliable  estimates  of  the  DKP 
and  Reynolds  stresses  depend  upon  the  accuracy  of  velocity 
estimates.  For  strong  currents  such  as  the  Gulf  Stream  and 
Kuroshio,  positions  should  be  determined  approximately  hourly 
with  a position  accuracy  of  the  order  of  450  meters. 


3 


4 


3)  The  position  fixing  system  should  be  designed  around 
minimizing  the  error  in  velocity  rather  than  acceleration. 

4)  The  Random  Access  Measurement  System  onboard  Nimbus  6 
represents  the  state  of  the  art  technology  for  position 
fixing  systems  and  data  telemetry.  It  is  adequate  for 
describing  scales  of  variability  of  the  order  of  24  hours, 
thus,  it  can  provide  reliable  estimates  of  velocity  and 
acceleration  but  not  for  DKP  and  Reynolds  stress. 

5)  Appendix  C breaks  down  the  problem  of  specifying  a posi- 
tion fixing  system  adequate  for  a strong  boundary  current 
kinematics  into  two  distinct  parts.  One  part  is  a specifi- 
cation of  the  position  fixing  system  based  on  positional 
accuracy,  rate  at  which  fixes  could  be  determined,  weight 
and  cost  per  buoy,  and  coverage  area.  The  following  sys- 
tems were  considered:  RAMS  onboard  TIROS  N,  OMEGA,  LORAN-C, 
the  Navy  Navigation  Satellite  System,  (HNSS)  and  the  Global 
Position  System  (GPS).  From  the  standpoint  of  position 
accuracy  alone,  the  CPS  appears  best  with  the  LORAN-C  next. 

6)  The  second  part  of  the  position  fixing  problem  considered 

in  Appendix  C was  the  telemetry  and  data  link.  The  following 
telemetry  systems  were  considered:  TIROS  N,  GOES,  MARISAT, 
ATS-5,  and  an  HF  telemetry  link.  Based  on  power  require- 
ments, cost,  range,  and  reliability,  it  appears  that  a 
TIROS  N telemetry  link  is  best. 


4 


B.  Recommendations 


1)  The  best  combination  of  position  fixing  system  and  telemetry 
link  for  a study  of  small-scale  variability  in  strong  bound- 
ary currents  appears  to  be  a LORAN-C  position  fixing  system 
with  a telemetry  link  through  the  TIROS  N.  The  cost  is 
expected  to  be  about  $6,000  per  unit.  This  incluses  hull 
costs  and  batteries  for  approximately  thirty  days.  The 
LORAN-C  was  picked  over  GPS  because  of  the  GPS  schedule  and 
uncertainties  in  user  costs,  hardware  costs,  and  power 
consumption. 

2)  Appendix  C briefly  describes  the  problem  of  drogue  slippage 
errors  of  drifters.  The  conclusion  was  that  there  were  no 
firm  data  which  reliably  correlates  the  buoy  motions  with 
the  effects  of  winds  and  waves.  It  is  recommended  that  the 
study  of  these  effects  be  initiated. 

3)  After  the  experience  with  Japanese  institutions  during  the 
pilot  experiment,  it  is  our  recommendation  that  it  is  cost 
effective  to  conduct  the  boundary  current  study  in  the  Kuro- 
shio  rather  than  the  Gulf  Stream.  The  following  reasons 
dictate  this: 

a)  The  Japanese  have  in  existence  an  extensive  monitoring 
program  of  the  Kuroshio.  They  seem  anxious  to  cooperate 
with  American  scientists  interested  in  studying  this 
current.  Thus,  it  would  be  possible  to  conduct  the 
experiment  with  Japanese  facilities  and  no  shiptime 
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cost  obligations  to  the  Office  of  Naval  Research. 


b)  The  extensive  monitoring  program  of  the  Kuroshio  has 
been  in  existence  for  many  years.  Such  an  historical 
data  base  does  not  exist  for  the  Gulf  Stream. 

c)  Once  the  feasibility  of  the  position  fixing  and  data 
relay  system  has  been  demonstrated  to  the  Japanese, 
it  is  expected  that  the  Japan  Hydrographic  Department 
would  underwrite  the  cost  of  the  drifters,  resulting 
in  further  savings  to  the  American  effort. 
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TRACKING  THE  KUROSHIO  BY  NIMBUS  6 


A.  C.  Vastano 


A.  D.  Klrwan,  Jr. 


Department  of  Oceanography 
Texas  A&M  University 
College  Station,  Texas  77843 


During  February,  1977,  we  had  the  opportunity  of  participating  with 
scientists  of  the  Japan  Hydrographic  Department  and  the  Ocean  Research 
Institute  of  the  University  of  Tokyo  in  a pilot  study  of  the  Kuroshio. 
This  is  a preliminary  report  of  our  efforts. 

As  one  of  the  major  Western  Boundary  Currents  (WBC) , the  Kuroshio 
has  long  been  an  object  of  intense  study.  In  many  respects,  the  Kuro- 
shio is  quite  similar  to  the  Gulf  Stream,  another  well-studied  WBC. 

There  is,  however,  one  major  difference.  The  Gulf  Stream  has  only  one 
* % 

preferred  path  as  it  flows  along  the  eastern  seaboard  of  the  U.S.,  but 
in  the  region  south  of  the  main  island  of  Honshu,  the  Kuroshio  follows 
one  of  two  possible  paths.  One  path  follows  quite  closely  the  coastline 
of  Honshu.  In  the  other,  a large  meander  develops  south  of  Honshu.  Thi3 
meander  is  believed  to  extend  as  far  south  as  30°N  or  about  300  km  from 
the  other  path.  Examples  of  these  two  paths  are  shown  in  figure  1 which 
were  taken  from  bulletins  that  are  distributed  regularly  by  the  Japan 
Hydrographic  Department. 

When  the  Kuroshio  is  in  its  meander  mode,  a cold  water  mass  develops 
between  the  meander  and  the  coastline.  South  of  the  meander,  a warm  core 
eddy  is  believed  to  exist.  This  meander  circulates  some  of  the  Kuroshio 
water  with  the  subtropical  counter  current  system  (Rikiishi  and  Yoshida, 
1974)  and  perhaps  the  Kuroshio  itself  south  of  Kyushu. 

Because  the  strong  currents  of  the  Kuroshio  have  such  a significant 
effect  on  the  shipping  industry,  and  the  intrusion  of  the  cold  water 
mass  greatly  affects  the  fishing  industry  and  weather  over  the  northern 
part  of  Honshu,  Japan  expends  a considerable  effort  in  monitoring  the 
Kuroshio.  This  includes  almost  continuous  ship  surveys  by  such  diverse 
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agencies  as  Japan  Fisheries  Agency,  Japan  Meteorological  Agency,  and 
the  Hydrographic  Department.  These  ship  surveys  include  closely  spaced 
XBT  stations,  some  hydrographic  stations,  and  continuous  underway  GEK 
measurements.  The  Hydrographic  Department  is  responsible  for  assimi- 
lating all  of  the  data  from  these  diverse  sources  and  producing  a bul- 
letin on  the  location  and  intensity  of  the  Kuroshio.  This  bulletin 
comes  out  every  other  Saturday  and  includes  observations  made  as  late 
as  the  previous  day.  Figuring  conservatively,  over  six-hundred  days  of 
ship  time  per  year  are  allocated  by  the  Japanese  government  for  moni- 
toring the  Kuroshio.  Despite  this  effort,  the  Japanese  feel  that  there 
are  significant  scales  of  motion  of  the  Kuroshio  which  are  not  detected 
by  either  the  frequency  or  the  scale  of  the  surveys. 

Because  one  of  us  had  previously  successfully  tracked  the  Gulf 
Stream  (Kirwan,  et  al. , 1976)  by  the  Nimbus  6 satellite,  the  Japanese 
were  anxious  to  test  this  technique  with  the  Kuroshio.  Thus,  we  were 
invited  to  participate  in  a cooperative  pilot  experiment  in  February, 
1977. 

We  deployed  four  drifters,  each  drogued  at  100  m by  a 9.2  m person- 
nel parachute,  off  the  island  of  Kyushu.  The  deployments  were  made  from 
the  R/V  TAKUYO  which  was  made  available  by  the  Hydrographic  Department 
of  the  Marine  Safety  Agency  of  Japan. 

Tills  pilot  experiment  had  a number  of  goals.  First,  we  wished  to 
compare  the  results  of  the  drifter  tracks  and  velocities  with  that 
obtained  by  GEK  and  the  hydrographic  surveys.  The  second  goal  was  to 
verify  the  existence  of  the  cold  and  warm  water  eddies.  The  third 
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objective  was  to  study  the  variability  of  the  currents  in  the  extension 


region. 

The  drifters  are  tracked  by  the  Random  Access  Measurement  System 
(RAMS)  onboard  Nimbus  6.  A brief  summary  of  the  way  the  system  operates 
is  as  follows:  The  transmit  terminal  in  the  buoy  transmits  approximately 
1 second  every  minute.  The  transmission  contains  a four-digit  identifier 
and  four  eight-digit  data  words.  Actually,  the  system  has  the  capability 
of  transmitting  eight  data  words,  four  words  in  alternate  transmissions. 
The  data  are  recorded  by  the  satellite  during  the  fly-by  (^20  min.)  . 

The  buoy  positions  are  determined  by  Doppler  shifts  of  the  transmissions 
occurring  during  the  fly-by. 

The  satellite  is  in  a sun-synchronous  polar  orbit  which  crosses 
the  equator  northward  bound  at  local  noon  and  southward  at  midnight. 

The  orbital  period  is  about  107  minutes.  Each  orbit  crosses  the  equator 
27°  longitude  further  west  than  the  previous  orbit. 

Discussion  of  Trajectories  West  oE  14fl°E 

The  four  drifters  were  deployed  along  a line  normal  to  the  axis  of 
the  Kuroshio  just  east  of  Honshu.  The  trajectories  up  to  Julian  day  110 
are  shown  in  figure  2.  Drifter  106  was  deployed  just  to  the  landward 
side  of  the  high  speed  core  of  the  Kuroshio  as  determined  by  GEK. 

Drifter  130  was  deployed  just  to  the  seaward  side  of  this  core.  Drifters 
307  and  341  respectively  were  deployed  out  on  the  flank  of  the  current. 
The  deployment  scheme  was  selected  with  the  expectation  that  106  would 
break  off  into  the  cold  core  eddy  and  that  341  would  peel  off  into  the 
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warm  core  eddy. 


Figure  2 shows  that  the  drifters  all  moved  at  approximately  the 
same  velocity  (2  - 2%  knots)  until  Julian  day  54.  At  that  time,  341 
took  off  to  the  southeast  along  the  meander  axis.  However,  as  indi- 


cated in  figure  2,  it  did  not  leave  the  Kuroshio  as  we  had  anticipated. 
At  the  other  extreme  was  106,  which  from  day  56  to  day  60  was  entrained 
in  a small  anticyclonic  eddy.  It  rejoined  the  main  flow  of  the  Kuro- 
shio on  day  62. 

This  is  a puzzling  result.  The  model  of  White  and  McCreary  (1976) 
suggests  that  106  should  have  undergone  the  greatest  acceleration  as 
the  flow  enters  the  meander  and  that  341  would  exhibit  the  least  acce- 
leration. Instead,  there  appears  to  be  an  area  of  stagnation  on  the 
northern  part  of  the  Kuroshio  as  it  starts  the  meander. 

Examination  of  figure  2 shows  that  the  trajectories  of  all  four 
drifters,  along  the  eastern  flank  of  the  meander,  are  quite  close, 
even  though  there  is  as  much  as  eleven  days  difference  in  the  time  of 
traverse.  This  supports  the  view  of  Japanese  scientists  that  during 
days  47  - 59  the  meander  was  in  an  extreme  western  position,  and  that 
it  was  beginning  to  move  eastward.  Thus,  in  eleven  days  the  western 
part  of  the  current  could  have  moved  20  km  to  the  east.  Such 
east-west  migrations  of  the  meander  are  well  documented  in  the  Japanese 
oceangraphic  literature. 

Figure  3 compares  the  four  trajectories  with  the  Japanese  analysis 
of  the  current  location  and  the  temperature  field  at  100  m for  days 
47  - 59.  The  dashed  lines  in  these  figures  are  areas  in  which  there  is 
a paucity  of  hydrographic  or  GEK  data.  The  agreement  with  the  current 
chart  generally  is  quite  good  with  only  minor  discrepancies  where  there 
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is  a lack  of  data.  Note  that  on  the  western  leg  of  the  meander,  two 


days  are  missing  in  the  track  of  307,  thus  the  straight  line  between 
date  points  is  not  a good  representation  of  the  actual  track.  Finally, 
there  is  no  evidence  in  the  hydrographic  or  current  charts  of  the 
anticyclonic  eddy  detected  by  106,  but  it  is  of  such  small  scale  that 
it  might  not  have  been  detected  by  the  survey. 

Data  from  days  60  through  74  are  compared  in  figure  4.  Up  to  the 
time  the  drifters  reach  the  Izu  Ridge,  which  rims  southward  along  140°W, 
the  100  m temperature  field,  the  current  chart,  and  trajectories  are  in 
excellent  agreement.  It  was  encouraging  to  note  that  307  went  straight 
through  the  Izu  Ridge  right  in  the  core  of  the  current.  The  track  of 
130  as  it  peels  off  from  the  main  flow  east  of  the  ridge  is  also  indi- 
cated in  the  current  chart.  Drifter  341,  after  it  reaches  the  ridge, 
enters  an  anticyclonic  eddy  which  is  not  inconsistent  with  the  surface 
flow  field  as  determined  by  GEK. 

For  days  75  -90,  only  106  and  341  remain  in  the  survey  area.  The 
latter  is  seen  in  figure  5 to  be  following  the  main  current  axis  east 
of  the  Izu  Ridge.  Drifter  341  drifted  slowly  southward  (about  .1  kt.) 
until  it  reached  30°  30'.  At  this  point  it  turned  eastward  and  then 
northward.  This  track  is  quite  consistent  with  the  GEfC  survey. 

Discussion  of  Trajectories  in  Kuroshio  Extension 
It  is  well  known  that  the  flow  east  of  the  Izu  Ridge  is  quite 
complex.  Data  from  numerous  hydrographic  surveys  indicate  a very 
complicated  structure  of  the  dynamic  topography  with  rapid  large-scale 
changes.  Ocean  fronts  and  large  eddies  are  conspicuous  features  here. 
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The  trajectories  in  this  region  (see  figure  2)  attest  to  this. 
Drifters  307  and  130  came  through  the  Izu  Ridge  at  the  same  point  but 
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separated  in  time  by  about  five  days.  Their  initial  tracks  upon 
departing  the  ridge  are  nearly  identical  to  latitude  148°.  At  this 
point,  however,  130  broke  off  in  a large  warm  core  eddy  north  of  the 
Kuroshio.  Drifter  307,  on  the  other  hand,  follows  a wave-like  pattern 
to  approximately  151°  longitude,  at  which  time  it  did  an  about-face 
and  started  flowing  to  the  southwest.  On  day  95  it  turned  abruptly 
northward  and  then  on  day  105  eastward. 

Drifters  106  and  341  have  been  moving  slowly  just  to  the  south  of 
the  mean  axis  of  the  Kuroshio.  They  arrived  in  this  area  by  completely 
different  routes;  341  having  performed  an  end  run  around  the  ridge, 
while  106  came  shooting  through  the  gap  at  33°  30’.  Since  day  109, 
however,  341  has  been  following  closely  the  trajectory  of  307,  except 
that  it  looped  much  further  to  the  north. 

A more  detailed  interpretation  of  these  tracks  cannot  be  made 
without  concurrent  hydrographic  or  infrared  satellite  data.  The  hydro- 
graphic  data  is  being  obtained  through  the  North  Pacific  Experiment 
(NORPAX)  TRANSPAC  XBT  program.  It  will  be  several  months,  however, 
before  it  is  available  to  make  such  a comparison.  It  should  also  be 
noted  that  data  from  the  drifters  is  still  being  received  from  N-i'ous  6. 

Summary 

To  summarize,  the  agreement  between  the  trajectories  and  the 
surveys  is  excellent.  Main  features  of  the  Kuroshio  such  as  path  and 
velocity  are  duplicated  in  both  the  drifter  and  survey  data  sets.  The 
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trajectories,  however,  failed  to  detect  either  the  warm  or  cold  core 
eddies  associated  with  the  meander.  Also,  it  was  quite  remarkable  to 
find  that  130  pinched  off  into  a warm  core  ring. 

From  the  Navy’s  standpoint,  one  of  the  most  significant  results 
is  the  demonstration  of  a (nearly)  real  time  data  return  capability 
provided  by  Nimbus  6.  The  RAMS  on  this  satellite  ha the  capability 
of  returning  not  only  position  but  environmental  data  as  well.  The 
transmit  terminal  electronics  are  very  reliable  and  draw  very  low  power 
so  that  long-term,  remote,  and  inexpensive  sensing  capability  is  possible 
now  for  operational  use  by  the  fleet. 
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Figure  Captions 


Figure  1 - Two  possible  paths  of  the  Kuroshio  as  shown  from  Japan 
Hydrographic  Department  bulletins:  (a)  1 June  - 17  June,  1975 
(b)  from  1 December  - 15  December  1975. 

Figure  2 - Summary  of  drifter  trajectories  from  time  of  deployment 
to  Julian  day  125  of  1977. 

Figure  3 - Comparison  of  drifter  tracks  with  current  (a)  and 

temperature  (b)  for  days  47-59  as  provided  in  Hydrographic  Department 
bulletin. 

Figure  4 - Comparison  of  drifter  tracks  with  current  (a)  and 

temperature  (b)  for  days  60-74  as  provided  in  Hydrographic  Department 
bulletin. 

Figure  5 - Comparison  of  drifter  tracks  with  current  (a)  and 

temperature  (b)  for  days  75-90  as  provided  in  Hydrographic  Department 
bulletin. 
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Coastline  path  of  the  Kuroshio  Current 


Meander  path  of  the  Kuroshio  Current 


Summary  of  temperature 


K I _2C  I ,0£  I 


Summary  of  surface  current  for  the  time  period  Julion  day  60  to  74 


temperature  (C°)  at  100  m for  the  time  period  Julian  day  60  to  74 


140 


Summary  of  surface  current  for  the  time  period  Julian  day  73  to  90 
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ABSTRACT 


This  is  a report  of  the  effect  of  random  position  error  and 
the  rate  of  position  fixing  on  the  inference  of  velocity  and 
acceleration  from  drifter  trajectory  data.  The  analysis  focuses 
on  the  relation  between  these  two  system  parameters  and  the 
natural  ocean  time  scales.  Direct  effects  of  ocean  space  scales, 
however,  are  ignored.  It  is  shown  that  for  a given  time  scale 
there  is  an  optimum  combination  of  position  fixing  sampling  rate 
and  position  error.  As  intuitively  expected,  the  limiting  factor 
for  sampling  rates  higher  than  the  optimum  is  position  error.  A 
more  serious  problem  arises  when  the  sampling  rate  is  less  than 
the  optimum.  Then  it  is  shown  that  the  limiting  effect  on  errors 
in  velocity  and  acceleration  is  aliasing  by  ocean  motions  with 
time  scales  less  than  the  sampling  rate. 

The  analysis  results  are  used  to  evaluate  four  position  fixing 
systems  previously  reported  in  the  literature,  as  well  as  the  ARCOS 
system,  which  will  be  in  use  during  the  First  Carp  Global  Experiment 
(FGGE)  . Tli is  shows  the  Over  the  Horizon  Radar  (0111!)  system  to  be 
the  least  effective  and  local  tracking  by  radar  to  be  the  most 
effective.  The  Random  Access  Measurement  System  (RAMS)  and  Eolc 
appear  to  be  reasonable  compromises  between  these  two  extremes; 
however,  there  are  some  ocean  time  scales  which  cannot  be  investi- 
gated by  these  systems.  The  ARGOS  system  will  be  a significant 
advance  over  these  last  two  systems. 
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1.  Introduction 


A recent  area  of  interest  in  oceanography  has  been  the  de- 
scription of  die  horizontal  scales  of  the  velocity  field  of  the 
ocean  through  observations  of  drifters.  For  example,  Freeland 
et  al.  (1975)  have  presented  some  intriguing  results  of  deep  flow 
in  the  western  North  Atlantic  from  Sofar  floats.  Earlier  work  in 
this  area  has  been  performed  by  Swallow  and  Worthington  (1961) 
and  Rossby  and  Webb  (1971).  Dickson  and  Baxter  (1972)  and 
Cresswell  (1976)  have  reported  some  current  observations  from 
satellite-tracked  drifters.  Kirwan  e_t  a_l.  (1976)  and  Richardson 
(1976)  also  have  employed  this  technique  Cor  studying  the  Gulf 
Stream  and  cyclonic  eddies  in  the  North  Atlantic.  Kirwan  and 
McNally  (1975)  used  Stanford  Research  Institute  Over  the  Horizon 
Radar  (UHR)  for  studying  the  North  Pacific  Cm  tent.  Molinari 
and  Kirwan  (1975)  determined  the  differential  kinematic  proper- 
ties (DKP)  in  the  vorticity  balance  in  the  Yucatan  Current  by 
the  use  of  radar-tracked  drifters.  A similar  technique  has  been 
used  by  Stevenson  ct  al.  (1974)  in  studying  upwelling  off  the 
coast  of  Oregon. 

The  studies  cited  above  have  determined  a number  of  kinematic 
properties  of  the  ocean.  These  include  trajectories,  velocities, 
Coriolis  and  inertial  accelerations  and  the  DKP  such  as  divergence 
and  vorticity.  The  sampling  rate  and  position  accuracy  of  the 
position  fixing  systems  employed  in  the  cited  studies  have  varied 
tremendously.  However,  an  analysis  of  how  these  characteristics 
affect  the  accuracy  of  the  calculations  has  not  yet  been  made. 
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Kirwan  et  al.  (19  76)  estimated  the  accuracy  of  estimates  of 
velocity  and  acceleration  from  the  Random  Access  Measurement  System 
(RAMS)  by  regarding  successive  position  measurements  as  a random 
time  series  in  which  the  position  uncertainty  <3  is  the  standard 
deviation  of  the  position  measurement.  Standard  deviations  for  ve- 
locity and  acceleration  wore  obtained  by  familiar  formulas  for  dif- 
ferencing random  variables.  Obviously,  such  an  approach  cannot  be 
applied  to  flow  in  which  there  may  be  significant  natural  and 
regular  variability  whose  time  scales  are  comparable  to  the  sampling 
period.  An  example  is  sinusoidal  flow  whose  period  is  less  than  the 
sampling  period. 

Okubo  and  Ebbesmeyer  (1976)  have  investigated  a number  of 
theoretical  problems  involved  in  estimating  kinematic  quantities 


from  drifter  data.  However,  the  effects  of  position  error  and 
sampling  rate  on  the  calculations  were  outside,  the  scope  of  their 
analysis . 

The  purpose  here  is  to  extend  these  last  two  studies  so  that 
effects  of  position  error  and  sampling  period  are  assessed 
properly.  It  is  hoped  that  the  analysis  will  result,  not  only  in 
a better  feel  for  the  accuracy  of  previous  calculations,  hut  also 
as  a guide  for  the  design  of  new  position  fixing  systems,  as  well 
as  criteria  for  planning  experiments  to  determine  the  horizontal 
scales  of  variability  in  the  ocean. 

It  is  recognized  that  there  are  other  factors  which  may  affect 
the  quality  of  drifter  data.  For  example,  with  surface  drifters. 
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winds  and  surface  currents  can  produce  considerable  errors  in  the 
trajectories  and,  perhaps,  the  velocities.  These  problems  are 
treated  elsewhere  (Kirwan,  ej^  al . , 1973). 

2.  Errors  in  velocity  and  acceleration 

The  basic  kinematic  information  obtained  from  drifters  is 
trajectories.  Other  kinematic  properties  such  as  velocity  and 
acceleration  are  inferred.  'Typically,  the  trajectory  data  is 
obtained  at  discrete  times  with  each  position  observation  subject 
to  a random  error.  A schematic  of  a typical  trajectory  with 
position  error  bars  is  shown  in  Figure  1.  The  observed  position 
at  t^  is  denoted  as  X^.  This  is  decomposed  into 


Xi  = 7i  + ri 

Here  is  the  expected  or  Liue  position,  and  r^  is  the  position 
error.  For  our  purposes  the  position  error  will  be  taken  as 
random  with  the  following  properties: 


rx  = 0, 


r . r . 
i J 


j2  6 


(2) 


ij 


Equation  (2)  asserts  that  the  random  position  error  has  zero  mean 
and  that  it  is  uncorrelated  with  the  position  error  at  other  times. 
However,  the  variance  of  the  position  error  at  any  one  time  is  the 
same  as  it  is  at  any  other  time.  Also,  for  convenience  it  will  be 
assumed  that  the  position  sampling  interval  is  constant  and  is 
denoted  by 


for  all  i. 
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The  true  position  at  times  t^+^  and  can  8iven  as  a 
Taylor's  Expansion  in  time  about  the  true  position  at  time  t^. 
This  is 


zi+l  - 

Ai 

1 o 1 

+ V.  A + 7A.  + yr  A. 

A3  + — 
4! 

A. 

l 

A4  +... 

(3) 

zi-l  + 

z. 

1 

- vi 4 + K 42  - h \ 

A3  + 

« « 

A. 

i 

A4  +. . . 

(4) 

In  these  equations  is  the  velocity  of  the  drifter  at  time  t^, 

and  A.  is  its  acceleration.  A.  and  A.  are  the  first  and  second 
r 1 i 

derivatives  of  the  acceleration.  They  are  termed  respectively 
jerk  and  quake.  All  of  these  quantities  are  evaluated  at  the 
true  position  (7..)  and  not  the  observed.  In  practice,  only  esti- 
mates of  these  quantities  at  the  observed  positions  (X^)  are 
available . 

First  consider  the  errors  associated  with  the  estimates  of 
velocity  and  accelerations  from  individual  drifters.  A velocity 
estimate  at  time  t^  is  obtained  from  the  centered  difference  formula 

xtn  - xi-i  ..  zm  - Vi.  + rt+i  - ri-i  (5 

2A  2A  * 

Equation  (5)  represents  the  average  of  the  velocity  in  the  interval 

h-i  to  fci+i. 


R-5 


Inserting  (3)  and  (4)  into  (5) , we  obtain 


X.  , - X.  , A.  r...,  - r,  , 

-~Va  ’ ' - Vi  + 3T  *2  + 2A  '--+°  (A4) 


(6) 


Squaring  and  averaging  this  last  result  gives  the  mean  square  ve- 
locity error: 


X.,,  - X.  . 

e ■ - v 

v l 2 A 


t a a4 

i / ~ (3!)^ 


+ a2/ 2A2  + 0 (A6) 


(7) 


The  acceleration  estimate  at  time  t_^  is  given  by 


X.M  - 2X.  + X 

1 1 -r-i = A,  + ~ A^  + 


. , .v.  r.,.  - 2r . + r.  . 

1 J"1  ■ ■ JL  A2  x -111 1 1=1  4.  nr a4 


A, 

^ “i  ■ 12  T A*" 

Squaring  and  averaging  this  last  result  yields 


I-  0 (A4)  (8) 


E = 
a 


2 

(M 

Aij 

(9) 


It  is  seen  from  (7)  and  (9)  that  for  small  sampling  periods  the 

-2 

position  error  term  dominates and  the  errors  are  proportional  to  A' 

-4 

and  A for  velocity  and  acceleration,  respectively.  This  is  the 
error  assumed  by  Kirwan  et_  aK  (1976).  On  t'ac  other  hand,  for  large 
sampling  periods,  the  position  error  term  is  small,  but  the  higher  order 
terms  in  (3)  and  (4)  dominate.  Thus,  for  large  A both  errors  go  as  A1*. 
In  both  cases  then  there  are  A's  which  give  a minimum  mean  square 
error.  Table  1 summarizes  the  optimum  A's,  denoted  A^,  and  the 


associated  errors. 


Table  1 


Statistics 


E 


v 


E 

a 


(3o/A.)1/3 

61/8[l2a/A.]1/4 


Minimum  Value  of  Error 


3p1a2/3.'J2/3  4 2/3 
aA . / v/  6 


These  formulae  show  that  A^  is  determined  solely  by  the  variance 
of  the  position  error,  which  is  a position  fixing  system  charac- 
teristic, and  the  higher  order  nonlinear  terms  in  the  Taylor's 
expansion,  which  is  an  ocean  characteristic.  The  table  also  shows 
that  the  minimum  rms  error  increases  with  the  increasing  rms  posi- 
tion error. 

Figure  2 is  a plot  of  a vs  A for  both  velocity  and  acceleration. 

In  this  graph  A is  varied  from  10  ' to  10  3 km  dy  3 and  A from  1 03  to 
~*2  “A 

10  km  dy  . (We  have  found  it  convenient  to  use  km  and  dy  as 
space  and  time  scales  for  large  scale  oceanographic  problems.  The 
conversion  to  these  units  from  m and  s for  velocity,  acceleration, 
jerk  (A),  and  quake  (A),  respectively,  are  1 ms  3 “ 86.4  km  dy 
1 ins-2  = 7.46-106  km  dy"2 , 1 ms"3  = 6. 45  *10 11  km  dy"3, 

1 ms  * = 5.57* 10 34  km  dy  ^ . ) 

The  Am's  arc  especially  signif icant:  in  that  they  represent  a 

natural  division  between  two  types  of  analyses  that  can  be  made  on 

the  observations.  For  A > A the  finite  difference  algorithms 

m 

given  above  are  efficient  techniques  for  estimating  velocity  and 
acceleration.  However,  for  A < Am  these  formulae  are  not  adequate. 
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Rather,  it  is  better  to  employ  specially  designed  filters  for 
calculating  velocities  and  accelerations  from  the  observations. 

Of  course,  knowledge  of  a2  and  other  statistical  properties  of 
the  position  error  can  be  used  in  the  construction  of  such  fil- 
ters. Such  techniques  can  reduce  the  error  to  less  than  that 
specified  by  the  A^. 

The  algorithms  yield  velocity  and  acceleration  estimates  as  a 
function  of  time.  From  the  Eulerian  standpoint,  the  nonrandom 
variability  in  these  estimates  is  the  result  of  both  spatial  and 
temporal  effects.  By  adhering  to  a purely  I.ograngian  description 
here  the  difficult  question  of  relating  the  effects  of  position 
error  and  sampling  rate  to  horizontal  variability  is  avoided. 

At  this  point  it  is  logical  to  question  the  effect  that  posi- 
tion error  and  sampling  period  have  on  the  calculations  of  the  1)KI* 
from  drifter  clusters.  Viewed  in  the  abstract  this  is  a formidable 
problem  because  the  equations  from  which  the  DKB  are  calculated 
involve  ratios  of  velocity  moments  with  respect  to  the  center  of 
mass  of  the  array  and  second  order  position  moments.  However,  with 
the  stipulation  that  the  estimated  center  of  mass  can  be  replaced  by 
the  true,  then  it  can  be  shown  that  the  D'tCP  are  merely  linear  func- 
tions of  the  velocities.  From  this  it  follows  (Cramer,  1066)  that 
the  mean  square  errors  (arising  from  position  uncertainty  and  sampling 
period)  in  the  DKP  are  proportional  to  the  mean  square  velocity  errors. 
The  appendix  gives  details  of  this  derivation  along  wiLh  a discussion 
of  the  ramifications  of  the  simplifying  assumption. 


3.  Assessment  of  existing  and  proposed  position  fixing  systems 
The  different  position  fixing  systems  reported  in  the  litera- 
ture, or  planned,  cover  a broad  range  of  position  error  and  sampling 
rate.  Table  2 is  a representative  sample. 

Table  2 

Position  error  and  sampling  rate  for  various  position  fixing  systems 
Svs  tem/Inves  tigator  g(kn)  A (day) 


Local  Radar/S teve.nson  ct_  a_l.  (1974) 

Molinari  and  Kirwan  (1975) 

-2 

5-10 

10 

Eole/D.ickson  and  Baxter  (1972) 

Cresswell  (1976) 

1 

• 5 

RAMS /Kirwan  _et_  al . 

1 

.5 

OIIR/Kirwan  and  McNally  (19  75) 

25 

3 

ARGOS/ 

1 

.1 

Eole,  RAMS , and  ARGOS  arc  all  polar  orbiting  satellite  systems. 
The  first  has  been  phased  out,  the  second  is  operational,  and  the 
last  is  scheduled  to  be  in  1978.  ARGOS  differs  in  two  significant 
ways  from  the  previous  systems.  First,  its  position  accuracy  will 
probably  be  less  than  1 km.  Secondly,  it  will  have  two  satellites 
in  orbit  rather  than  one.  This  will  greatly  increase  the  number  of 
fixes  available  from  this  system.  Moreover,  these  fixes  will  be 
more  uniformly  spaced  in  time  than  in  the  previous  systems. 

Since  the  satellite  systems  are  polar  orbiting,  the  number  of 
fixes  vary  tremendously  with  latitude.  The  numbers  given  in 
table  2 are  conservative  estimates  for  a mid-latitude  region. 

In  order  to  assess  the  capability  of  these  systems,  it  is  neces- 
sary to  have  estimates  of  the  jerk  and  quake  of  ocean  motions. 
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Unfortunately,  calculations  of  these  higher  order  derivatives  from 
current  meter  or  drifter  data  may  be  contaminated  by  measurement 
errors  and,  perhaps,  insufficient  sampling  rates. 

This  obstacle  can,  in  part,  be  overcome  by  prescribing  the  flow 
by  a simple  model.  Such  an  approach  is  desirable  in  that  extreme 
situations  are  easily  investigated.  To  this  end,  consider  a periodic 
velocity  field  specified  by  a frequency  m and  scale  velocity  U.  The 
jerk  and  quake  are  readily  found  to  be  of  the  order  of 

A ~ 

A ~ Uu3 

Two  cases  are  of  interest.  For  the  "fast"  ocean  case,  take  U 
as  8.64*10  km  dy  ^ and  m as  2tt  rad  dy  This  might  correspond  to 
a daily  cycle  in  a strong  current  such  as  the  Gulf  Stream.  At  the 
other  extreme  is  the  "slow"  ocean  as  typified  by  a mcsoscale  eddy 
whose  scale  velocity  and  period  of  rotation  are  8.64  km  dy  ^ and  20  dy, 
respectively. 

Table  3 summarizes  the  velocity  and  acceleration  A 's  for  these 
two  cases  for  each  of  the  position  fixing  systems  mentioned  above. 
Inspection  of  this  table  shows  that  for  both  the  fast  and  slow  ocean 
cases,  the  velocity  A will  always  be  less  than  the  acceleration  A . 
This  means  that  if  the  position  fixing  system  A is  close  to  the 
velocity  A , the  acceleration  is  "over  sampled".  This  is  a desirable 
situation  because  special  filtering  techniques  can  be  used  to  provide 
better  estimates  of  the  acceleration  than  would  be  available  if  the 
system  A was  the  acceleration  A^.  Note  that  this  is  a system  charac- 
teristic which  applies  only  when  the  A is  close  to  the  Am-  If  the  A’s 
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Table  3 Velocity  and  Acceleration 


are  less  than  the  A^'s  and  the  velocity  spectrum  decreases  with 
frequency,  then  the  acceleration  spectrum  will  be  whiter  than  the 
velocity  spectrum.  For  comparable  accuracy  this  would  require  an 
acceleration  A to  be  less  than  the  velocity  A. 

Figure  3 is  a plot  of  ratio  of  the  rms  velocity  error  to  the 
Velocity  scale  as  a function  of  A.  This  figure  was  constructed  by 
specifying  the  jerk  and  quake  from  the  fast  and  slow  ocean  cases. 
This  shows  the  two  error  regimes  mentioned  above.  The  solidus , 
with  a slope  of  4,  to  the  right  of  the  minimum  is  the  region  in 
which  the  ocean  variability  effects  (jerk  and  quake)  are  the 
dominant  error  sources.  If  the  position  fixing  system  A intersects 
this  portion  of  the  curve,  there  is  little  that  can  be  done  to 
reduce  the  error. 

Such  is  not  the  case  if  the  system  A intersects  the  dashed 
lines  whose  slope  is  -2  which  is  to  the  left  of  the  minimum. 

This  is  the  region  in  which  the  error  indicated  in  the  figure  can 
be  significantly  reduced  by  filtering  the  data. 

It  is  stressed  that  the  error  ratio  given  in  figure  3 is  the 
worst  case  estimate  in  that  the  true  velocity  was  taken  as  sinu- 
soidal. In  reality,  there  is  frequently  a significant  steady  or 
nonsinusoidal  component.  Incorporation  of  a steady  component  in 
the  true  velocity  will  reduce  the  error  accordingly. 

The  calculations  described  above  suggest  that  RAMS  is  not 
suitable  for  investigating  diurnal  oscillations  of  the  Gulf  Stream. 
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However,  it  could  certainly  be  used  to  study  longer  period  phenomena 
there.  Note  that  figure  3 shows  that  the  A for  this  system  makes  it 
ideal  for  investigating  mesoscale  eddies  or  longer  period  phenomena. 

Neither  RAMS,  Eole  nor  the  OHR  is  likely  to  be  adequate  for  cal- 
culating DKP  unless  a large  number  of  drifters  are  used.  To  achieve 
a reasonable  accuracy  it  may  be  necessary  to  use  A's  significantly 
less  than  the  A^’s.  A position  fixing  system  A of  the  order  of  .1 
a day  and  a of  10  ^ km  might  suffice  for  this.  The  ARGOS  system 
approaches  this. 
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Appendix 

Errors  in  Differential  Kinematic  Properties 
In  order  to  analyze  the  effect  of  position  error  and  sampling 
rate  on  the  DKP,  it  is  necessary  to  indicate  both  the  X and  Y 
components  of  position  as  well  as  the  time  and  drifter  number.  To 
this  end  the  coordinates  relative  to  the  center  of  mass  of  the 
cluster  for  the  m^  drifter  are  denoted  as 


X = Z + r i = 1,2  ; m=l,...n  (A-l) 

lm  1m  im 

Here  the  first  subscript  refers  to  the  coordinate  direction  (X  or  Y) 
and  the  second  to  the  drifter  number.  As  before,  the  true  position 
is  Z and  r is  the  random  position  error. 

Following  Molinari  and  Kirwan  (1975)  and  Okubo  and  Ebbesmeyer 
(1976),  the  velocity  components  of  the  m*"'1  drifter  of  a cluster  of  n 
drifters  can  be  expressed  as  a Taylor's  Expansion  (in  space)  about 
the  cluster  center  of  mass.  Thus 


U = U + X.  + 
m o dX^  xm 


3V 

V = V + X,  + 
m o im 


(m  = l...n;  sum  on  i = 1,2)  (A-2) 


In  (A-2)  U and  V can  be  calculated  by  the  algorithms  given  above 
m m 

3U  3V 

while  the  coefficients  U , V , tt;  and  — can  be  determined  by 

O o dA,  dX. 

i l 

least  squares  provided  there  are  three  or  more  drifters  (Molinari 

and  Kirwan,  1975;  Okubo  and  Ebbesmeyer,  1976). 

The  aim  here  is  to  determine  how  the  position  error  and  sampling 

rate  affect  the  accuracy  of  the  estimates  of  U , 3U/3X. , etc.  For 

J o i 
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this  it  is  only  necessary  to  consider  the  U component  equation. 

The  assumption  is  now  made  that  the  expansion  given  by  (A-2) 
is  valid  relative  to  either  the  observed  or  true  center  of  mass 
if  the  true  coefficients  are  used.  For  many  applications  this  is 
not  a severe  restriction.  The  variance  of  the  difference  between 
these  two  position  is  a2/ n.  Thus,  the  estimated  center  of  mass 
will  be  within  one  standard  deviation  of  the  true.  Moreover,  it 
can  be  made  as  close  to  the  true  as  desired  by  increasing  the 
number  of  drifters.  One  a-priori  condition  on  how  close  the  esti- 
mated to  the  true  should  be  is  that  it  should  be  less  than  the  rms 
relative  displacements  of  the  drifters  from  the  estimated. 

With  this  assumption  it  is  seen  that  only  the  left  hand  side 
of  (A-2)  need  now  be  considered  as  random.  Thus,  the  problem  is 
reduced  to  that  of  multiple  regression  with  nonrandom  independent 
variables  (see  Cramer,  1966). 


Hie  following  nonrandom  functions  are  introduced: 


• - T (x.  - x.)  (X.  - X.)  i,  j = 1,5 

ij  1 Jm  J 


x.  = - = y x. 

in  im 

m~l 

P - det  U.J 

Mj,  j = iC\  jth  cofactor  of 
along  with  the  random  vector 

' 2 ® (xi„  ' h> 

m=l 


(A-3) 

(A-A) 

(A-5) 

(A-6) 


Here  U is  the  expected  value  of  the  X component  of  the  center  of 


mass  velocity. 
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The  use  of  (A- 3)  through  (A-6)  allows  us  to  express  the  least 

. , „ DU  , DU 

squares  estimates  of  Uq,  — , and  — as 


n 

X 

Um 

(A-7) 

nH=l 

{si 

NlI  + S2 

U,2}/P 

(A-8) 

{si 

«21  +s2 

«22I/P 

(A-9) 

JilL- 

3X, 


ILL 

3X2 


Note  from  (A-7)  through  (A-9)  that  each  of  these  coefficients 
is  a linear  function  of  the  drifter  velocities  which  are  random. 
In  such  cases  it  is  readily  shown  that  the  expected  values  of  the 
estimates  are  merely  the  same  linear  functions  of  the  expected 
values  of  the  velocities.  From  ( A-2 ) it  is  seen  that  the  latter 
involves  the  true  values  of  the  coefficients.  By  inserting 
(A-3)-(A-6)  into  (A-8)  and  (A-9) , it  is  found  that 


in 

3X, 


in 

3X, 


U = U 
o o 


(A- 10) 


i w“i 

A second  aspect  of  the  linear  relation  is  that  the  covariance 
matrix  of  the  estimates  are  given  by 

(A-ll) 


(U  - U )2  = e2/n 
o o 


/ ~ \ 

/ ^ 

3IJ  DU 

/ su 

DU 

3X.  3X. 

3X. 

3X 

1 7 

l J 

7 

= c2Lij/nL,  (i , j - 1,2) 


Here  the  maximum  likelihood  estimate  of  e2  is 
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(A- 12) 


i n jntt  *1 2 

- - T.,  U - U - ~ X.  (sum  on  i = 1,2)  (A-13) 

n m=l  |_  m o DX^  im  ' 


iv 


Eq.  (A-10)  says  that  the  expected  value  of  the  estimates  of  the 
center  of  mass  velocity  and  the  velocity  gradients  are  the  true 
values.  Eqs.  (A-ll)  and  (A-12)  show  that  the  variance  of  the  center 
of  mass  velocity  and  the  velocity  gradient  covariance  are  propor- 
tional to  the  regression  variance.  From  (A-13)  this  is  seen  to  be 
composed  of  the  moan  square  velocity  error  (velocity  variance)  and 
quadratics  involving  the  estimated  velocity  gradients.  The  last 
set  of  terms  may  be.  assumed  to  be  independent  of  the  sampling  rate 
(Molinari  and  Kirwan,  1975).  Thus,  the  covariances  are  all  propor- 
tional to  the  mean  square  velocity  error,  and  their  dependency  on  A 
is  the  same,  except  for  a displacement  along  the  ordinate. 
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Figure  Captions 

Figure  1 Schematic  of  drifter  trajectory  data  obtained  at 
discrete  times. 

Figure  2 Plot  position  error  versus  velocity  and  acceleration  A^'s. 

* **  «.  1 

Jerk  (A)  and  quake  (A)  are  varied  from  10  to  10  km  dy 

2 5 — /j 

and  10  to  10  km  dy  , respectively. 

Figure  3 Plot  of  normalized  rms  velocity  error  as  a function 

-2 

of  A.  Position  error  a is  varied  from  10  km  to  10  km. 
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INTRODUCTION 


At  the  present  time  the  Euoy  Transmit  Terminal  (BIT)  of  the  NIMBUS-6 
satellite  has  demonstrated  the  simplicity,  reliability  and  cost  effectiveness 
of  deriving  oceanographic  data  from  satellite- tracked  drifters  (Kirwan,  1977 
and  Richardson  et  al,  1977).  Errployed  in  an  expendable  mode,  these  buoys  have 
measured  useful  data  for  periods  of  over  one  year  in  such  places  as  the  North 
Pacific  at  latitudes  to  45JN.  The  satellite  positioning  system  has  provided 
a position  fix  and  enabled  the  transmission  of  on  board  sensor  data  at  fre- 
quencies of  2 to  8 times  per  day,  with  2 transmissions  assured  when  the 
satellite  is  nearly  overhead  twice  per  day.  The  positional  accuracies  reported 
from  controlled  tests  are  approximately  * Km,  901  of  the  time  when  the 
satellite  is  overhead  twice  per  day  and  less  when  passing  near  the  horizon. 

In  pieinning  an  experiment,  an  investigator  can  really  only  plan  on 
position  data  every  12  hours.  At  this  rate  there  is  a risk  that  energy  at 
the  diurnal  or  shorter  periods  my  be  aliased  (i.e.,  Nyquist  frequency  of 
1 cycle/day  = l/2At,  At  = sample  period).  As  a result,  the  existing  NIMBUS-6 
(to  be  succeeded  by  the  TIROS-N  satellite  in  1978)  satellite  is  useful  in  the 
study  of  slcwly  changing  phenomenon  and  large  scale  circulation. 

Dr.  A.  D.  Kirwan  of  Texas  A.&  M.  University  (TAMU)  has  chosen  to  measure 
the  notions  of  western  boundary  currents  in  general  and  the  Gulf  Stream  and 
Kuroshio  currents  in  particular.  Because  of  the  extent  of  such  currents  and 
cost  considerations  the  method  of  investigation  will  be  by  remotely-positioned, 
expendable,  drogued  drifting  buoys.  It  has,  however,  been  found  by  Robinson  et  al 
(1974)  and  others  that  one  of  the  dominant  energy  modes  for  Gulf  Stream" 
meanders  is  at  approximately  the  diurnal  period.  The  buoy  transmit  terminal 
(BIT)  of  either  the  NIKEUS-6  or  TIROS-N  satellite  is  therefore  an  inappro- 
priate method  of  directly  deriving  sufficient  buoy  trajectory  data.  The 
classical  manner  of  handling  aliasing  in  dynamic  systems  is  to  sharply  attenuate 
energy  at  frequencies  above  tdie  Nyquist  frequency  by  electronic  or  mechanical 
filters.  Drifting  buoy  motion  cannot  be  filtered  in  a like  manner,  nor  would 
it  be  desirable.  This  is  the  very  information  being  sought. 

Alternative  schemes  for  economically  measuring  buoy  position  at  a 
higher  rate  than  a BIT  will  be  explored  in  this  writing.  In  addition,  the 


trade-offs  of  buoy  positional  accuracy,  frequency  of  position,  and  drogue 
current  locking  accuracy  (i.e. , slippage)  will  be  discussed. 

2.0  OCCANOGRAPIUC  AND  SYSTEM  DESIGN  OBJECTIVES 

The  highest  frequency  energy  reported  for  Gulf  Stream  Motions 
(Robinson  et  al,  1974)  was  at  a period  of  approximately  24  hours.  They 
also  indicate  that  there  is  an  appearance  of  energetic  small  scale  structure 
(period  ^2-4  days)  which  ray  possibly  disappear  into  large  scale  structures 
(period  vl  day) . Because  their  detailed  data  were  limited  to  a 180-km  swath 
along  the  Gulf  Stream  axis,  it  was  not  possible  to  trace  the  smaller  scale 
disturbances  over  a wide  enough  area  to  document  behavior  conclusively . There 
is  even  very  limited  evidence  in  their  data  of  structure  at  periods  of  less 
than  one  day. 

Sampling  theory  indicates  that  position  data  should  be  recorded  at  least 
twice  per  sliortest  period  in  order  to  record  unaliascd  trajectory  data  at  this 
frequency.  In  general,  though,  it  is  desirable  to  have  at  least  5-10  data  points 
for  the  shortest  period.  Therefore,  a nominal  design  goal  will  be  to  acquire  a 
position  data  point  every  hour.  For  econonr/  of  onboard  buoy  power  it  will  be 
initially  assumed  that  the  *3  c error  footprint  of  a buoy  position  fix  does  not 
overlao  the  footprint  of  a subsequent  fix.  If  it  was  later  found  that  there 
was  an  overabundance  of  power  and  data  capability,  fixes  could  be  taken  more 
frequently  and  averaging  done.  With  the  above  assumption,  the  maximum  tolerable 
positioning  error  for  useful  measurentants  is  given  by  the  relation: 


where  V is  the  minimum  expected  velocity,  f the  sampling  frequency  (assumed 
at  1 hr-1) , and  e the  3 o bound  of  an  error  footprint  in  the  direction  of  buoy 
travel.  For  western  boundary  currents  such  as  the  Gulf  Stream  or  the  Kuroshio 
the  near-surface  velocity  is  generally  cjreatcr  than  0.5  knots  (i.e.,  .9  km/hr) 
which  dictates  a positional  accuracy  cf  *450  m for  hourly  position  fixes. 

The  planned  lifetime  of  the  system  electronic  design  is  30  days.  This 
value  is  based  on  an  estimate  of  the  longest  period  of  time  that  a buoy  would 
remain  witliin  the  current  regime  of  interest.  All  power  budget  figures  which 
follow  can  be  linearly  scaled  about  this  nominal  lifetime  while  confidence  on 
the  acquisition  of  certain  radio  signals  varies  wTith  buoy  location  with  re- 
spect to  stations. 

It  is  poor  jud ament  to  place  a sophisticated,  precise,  radio  positioning 
system  aboard  a drifting  buoy  for  the  purpose  of  monitorinq  currents  when 
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the  buoy  itself  my  have  drogue  slippage  between  position  fixes  which  is 
in  excess  of  the  positional  accuracy.  In  such  cases  the  derived  trajectory 
data  can  be  meaningless  in  light  of  the  defined  goals.  On  the  other  hand, 
it  is  recognized  tint  even  for  the  best  design,  under  certain  extreme  weather 
conditions,  highly  erroneous  data  can  arise  and  siould  be  discounted.  This 
writing  will  attenpt  to  bound  the  slippage  error  as  a function  of  the  partic- 
ular buoy  design  and  environment  by  applying  existing  data  and  analytic 
approaches  to  estimating  buoy  forces  and  drogue  slippage.  Lastly,  suggestions 
will  be  made  as  to  required  parameters  to  be  measured  and  means  by  which 
first  order  corrections  can  be  made  for  drogue  slippage. 

3.0  BUOY  POSITIONING  SYSTEMS 

The  problem  of  automatically  determining  the  position  of  a drifting 
buoy  breaks  down  into  two  parts  - the  position  determination  and  the  means  of 
reporting  the  position.  This  section  will  address  the  various  means  of  posi- 
tion determination  (both  present  and  future) , irony  of  which  also  provide  a 
means  of  automatically  reporting  the  data.  The  next  section  will  cover  the 
existing  schemes  for  the  telemetry  or  otherwise  reporting  buoy  position. 

3.1  TIROS-N  (or  NIMBUS-6) '/  RAMS  Satellite  System 

As  mentioned  in  the  introduction,  the  Random  Access  Measurement  System 
aboard  either  the  NIMBUS-6  satellite  (presently  orbitting  but  working  only 
half-time  every  other  day)  or  the  TIROS-N  satellite  (mid  1978  launch  date) 
would  not  be  an  appropriate  method  of  deriving  buoy  position  data  because 
of  the  infrequent  reliable  fixes.  Furthermore,  section  2 also  shows  that  its 
positional  accuracy  does  not  meet  the  design  requirements.  The  RAMS  system 
could,  however,  provide  a means  of  resolving  positional  ambiguities  or  re- 
ceiver error  problems  in  a DURAN  C positioning  system.  For  example,  occasion- 
ally a LORAN  C receiver,  which  normally  tracks  the  third  cycle  crossing  in  a 
pulse,  may  erroneously  track  the  second  cycle  crossing,  leading  to  an  error  of 
10  microseconds  or  approximately  a 3 km  error.  If  the  RAMS  system  wore  used 
as  the  means  for  DORAN  C data  tranmittal  such  errors  oould  be  revealed.  Other- 
wise the  RAMS  system  by  itself  is  not  useful  as  a prime  buoy  positioning  system. 
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3.2  OMEGA  Radio  Positioning 

OMEGA  is  a very  long  range,  hyperbolic  radio  positioning  system  that 
works  in  the  Very  Low  Frequency  (VLF)  range.  Lines  of  position  are  generated 
by  phase-difference  measurements  of  a conposite  threc-f requency , time- 
shared  system  from  each  of  eight  stations  giving  worldwide  coverage.  Wien 
fully  operable  the  eight  stations  will  be  at  the  following  locations  which 
exhibit  up  to  8,000  to  10,000  Im  baseline  distances: 

(1)  North  Dakota 

(2)  Reunion  Island,  Indian  Ocean 

(3)  Trinidad 

(4)  Norway 

(5)  Japan 

(6)  Australia 

(7 ) Argentina 

(8)  Hawaii 

At  present  all  but  the  Australian  station  arc  operable. 

Each  station  transmits  cadi  of  the  following  frequencies  for  approxi- 
rrately  one  second  (continuous -wave)  every  ten  seconds  in  a prescribed  format 
at  a 10  lew  power  level: 

fx  = 10.200  kHz 

f2  = 11.333  kHz 

f3  = 13.600  kHz 

The  particular  timing  format  of  a given  transmitter  allows  its  unique  identi- 
fication. The  time  synchronization  between  stations  is  maintained  by  atomic 
clocks. 

The  conposite  or  three-frequency  format  is  helpful  for  two  reasons: 

(1)  increased  lane  widths 

(2)  increased  knowledge  of  signal  phase  propagation 

As  stated  earlier,  a corplete  OMEGA  receiver  monitors  the  phases  of  three 
pulsed,  continuous  waves  from  each  transmitter.  If  a single  frequency 
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receiver  were  employed  it  v.ould  find  that  the  signal  phase  relationship 
along  a station  baseline  would  be  tlie  same  every'  X/2  (X  = wavelength) . For 
the  10.2  kHz  signal  alone  this  would  occur  every  15  km  (i.e. , the  single 
frequency  lane  v.'idth) . Therefore,  if  a sudden  ionospheric  disturbance  (SID) 
should  occur,  arising  in  a matter  of  minutes  from,  a solar  flare,  the  accom- 
panying ionospheric  lowering  could  result  in  as  much  as  a 40  centicycle 
(i.e.,  40  CEC  = .4  cycle)  change  in  signal  phase  for  _,n  all-daylight  path 
(Research  Triangle  Institute, 1973) . Furthermore,  as  the  solar  terminator 
passes  across  the  path  between  transmitter  and  receiver  as  much  as  a 25  CEC 
change  can  occur  at  the  receiver  over  a period  of  30  minutes  (Beukers,  1973) . 
Such  rapid  Occurrences  might  easily  lead  to  loss  of  lane  count  if  constant 
signal  monitoring  were  not  maintained. 

In  an  effort  to  artificallv  widen  the  lanes,  in  order  to  lessen  the 
amount  of  position  monitoring  required,  three  frequencies  wore  introduced. 

In  essence  the  larger  lane  widths  are  new  determined  by  the  difference  or 
beat  frequencies  between  the  transmitter  frequencies.  Therefore,  the  dif- 
ference frequency  between  the  10.200  and  11.333  kHz  signals  dictates  a lane 
width  of  133  kin. 

Tine  second  reason  for  using  the  oonposite  OMEGA  system  is  that  by 
taking  advantage  of  the  available  phase  information  at  more  than  one  fre- 
quency it  is  possible  to  enhance  computerized  predictions  of  signal  pliase 
propagation.  Such  a potential  arises  from  the  fact  that  phase  perturtetions 
at  VLF  are  to  some  extent  frequency- dependant  (Pierce,  1972). 

Without  any  corrections  at  all,  the  positional  accuracy  of  an  OMEGA, 
receiver  can  be  on  the  order  of  10  to  15  km.  By  estimating  errors  intro- 
duced by  diurnal  variations  such  as  variable  propagation  velocities  as  a 
function  of  solar  angle  and  variable  ground  conductivity  it  is  possible  to 
increase  the  accuracy  to  better  than  *2  km  (daytime)  and  better  than  - 3.5  km 
(nighttime)  (Research  Triangle  Institute,  1973).  This  accuracy  is  not  good 
enough  for  the  drifting  buoy  program  envisioned  nor  does  it  perform  any  tetter 
than  the  positional  accuracy  of  the  BIT/RAMS  package  which  is  on  the  NIMBUS-6 
or  TIROS-N  satellite  system. 

The  accuracy  of  the  OMEGA  system  only  approaches  an  acceptable  level 
when  it  is  used  in  the  differential  mode.  That  is,  a nonitor  station  should 


5 


1 


be  irDunted  at  a fixed  location  and  its  recorded  positional  changes  sub- 
tracted from  that  of  the  buoy.  For  buoy-monitor  ranges  (i.e.,  correlation 
distances)  of  less  than  20  km  the  positional  accuracy  can  be  as  good  as 
± 450  m (i.e.,  3 CCC)  for  data  which  are  instantaneoulsy  updated.  If,  how- 
ever, the  correlation  distance  exceeds  400  km,  positional  errors  approach- 
ing 4.5  kn  (i.e.,  30  CDC)  can  be  expected  (Beukers,  1973).  Because  buoys 
tracking  western  boundary  currents  can  be  expected  to  be  on  the  order  of 
hundreds  of  kilometers  from  land  nesses  it  does  not  seem  possible  to  employ 
differential  OMEGA..  On  the  other  hand  it  does  not  even  seem  economically 
feasible  to  think  of  employing  monitoring  stations  unless  they  were  already 
installed  for  other  reasons.  It  will , however,  be  mentioned  later  that  it 
is  hoped  that  the  OMEGA  system  might  fit  into  a multi-purpose  modular  drifting 
buoy  positioning  system  that  could  be  taken  anywhere  in  the  world. 


3.3  LOPnAN-C  Radio  Navigation 

LORAN-C  is  a high-accuracy  radio  navigation  system  which  appears  to 
be  the  best  choice  for  tracking  drifting  buoys  in  the  Gulf  Stream.  The 
system  provides  accuracy  of  greater  than  500  m at  a range  of  ™00  Tb 
system  is  also  highly  reliable  and  there  is  a large  amount  of  commercial 
receiving  equipment  available  * 


3.3.1  Operation  of  LORAN-C  Chains 
Lines  of  Position 

A position  fix  is  obtained  using  LORAN-C  by  determining  the  time 
difference  between  the  arrival  of  a pulse  from  a master  station  .and  each  of 
two  slave  stations.  Each  time  difference  defines  a hyperbolic  line  of  posi- 
tion (LOP).  The  intersection  of  two  LOPs  corresponds  to  the  position. 

Group  Repetition  Interval 

All  transmitters  operate  at  100  kHz  and  transmit  groups  of  pulses 
at  an  interval  known  as  the  Group  Repetition  Interval  (GRI) . The  GRI  is  ex- 
pressed in  tens  of  microseconds  and  varies  from  4000  to  9999.  Each  chain  (a 
master  and  up  to  4 slaves)  has  its  own  GRI. 
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Coding  Delay 

Slaves  are  positioned  such  tliat  at  least  too  can  be  received 
throughout  the  desired  coverage  area.  'Do  prevent  ambiguity  between  slaves, 
each  slave  transmits  after  the  master  by  a time  interval  known  as  the  coding 
delay.  Coding  delays  are  chosen  such  that  throughout  the  coverage  area,  pulses 
will  always  be  received  in  the  same  order,  that  is  Master,  Slave  x,  Slave  y. 

Slave  z,  etc.  This  is  s’nov.n  in  Figure  1. 

Blink  Code 

Each  group  of  pulses  consists  of  8 pulses,  with  the  exception  of  the 
Master,  which  transmits  9 pulses.  The  extra  pulse  serves  to  identify  the 
Master.  Additionally,  this  ninth  pulse  in  the  Master  group  may  be  switched 
on  and  off  bo  inform  users  of  a transmitter  malfunction.  Different  on-off 
patterns  are  used  to  tell  which  transmitter  is  malfunctioning.  This  is  loicwn 
as  the  blink  code. 

3.3.2  LORAN-C  Accuracy  and  Repea  debility 

High  accuracy  of  LOPAN-C  is  obtained  by  taking  the  time  of  arrival  to 
be  the  zero-crossing  at  the  end  of  the  third  cycle  of  the  groundwave , as 
shown  in  Figure  2.  The  shape  of  the  pulse  allows  this  third  cycle  zero-cross- 
ing to  be  identified.  Groundwave  propagation  is  stable  with  respect  to  time 
of  day,  year,  sunspot  cycle,  aid  ionospheric  conditions  in  general.  Skvwaves, 
reflected  off  different  levels  of  the  ionosphere,  will  arrive  between  35  usee 
and  1000  usee  after  the  corresponding  groundwave.  The  earliest  case  will  not 
interfere,  as  the  msasureiront  is  made  30  usee  after  the  beginning  of  the  pulse. 
Therefore,  the  groundwave  signal  can't  be  contaminated  by  its  own  sky  wave. 

However,  a pulse  may  be  contaminated  by  the  skywave  of  the  previous  pulse  if 
the  skywave  delay  approaches  1000  usee.  Automatic  receivers  can  identify 
these  delayed  pulses,  since  the  phase  of  the  100  kHz  carrier  is  changed  in 
each  pulse  in  a group  according  to  a prearranged  sequence  known  as  the  phase  code. 

The  accuracy  of  the  LORAh'-C  system  is  a function  of  tie  receiver  loca- 
tion witlin  the  coverage  area  as  well  as  time  measurement  noise, 
earth  conductivity  changes,  atmospheric  refraction  (for  cases  of  skywave  inter- 
ference) , and  system  synclironization . The  accuracy  as  a function  of  location, 
referred  to  as  cccrmtric  dilution  of  position  (GX)P) , can  become  rather  large 
for  hyperbolic  systems,  especially  in  regions  of  transmitter  baseline-extension. 
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Figure  1 - Exairple  of  Received  LDRAN-C  Signal 

(from  LORAN-C  User  Handbook) 


Figure  2 - Typical  LORAN-C  Pulse 

(from  LORAN-C  User  Handbook) 
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In  general,  the  overall  system  accuracy  is  within  500  meters,  95%  of  the  time; 
with  a relative  or  repeatability  error  of  15  to  100  meters  - depending  on  the 
above  conditions  (Fried,  1974).  In  practice,  the  measured  time  differences  at 
the  same  spot  will  vary  by  about  O.lusec.  Tiiis  time  difference  variation  re- 
sults in  positional  variations  which  depend  on  the  location  of  the  rcc'rver, 
master,  and  slaves.  Figure  3 shows  this  positional  variation  in  meters  for  the 
region  of  the  United  States  East  Coast  drain  (Micrologic  Inc. , 1975) . 

3.3.3  LORAN-C  Coverage 

The  existing  U.  S.  East  Coast  Chain  provides  adequate  coverage  for 
a Gulf  Stream  monitoring  experiment  out  to  a position  north  and  east  of  Ber- 
muda. Additionally,  coverage  is  available  in  many  other  parts  of  the  world, 
as  shown  in  Figure  4.  All  drains  transmit  on  the  same  frequency  (100  kHz) , 
but  generally  should  not  interfere  with  eadi  other,  as  all  drains  have  different 
GRIs.  Therefore,  the  same  receiving  equipment  co-old  be  used  with  any  existing 
chain. 

When  the  Coast  Guard  has  completed  installation  of  the  newest  stations, 
new  IORAN-C  drains  will  be  created,  leading  to  increased  coverage  as  shewn  in 
Figure  4.  This  dunce  will  mean  that  in  a few  locutions  such  as  the  southern 
part  of  tire  U.  S.  cast  coast  there  are  potential  interference  problems  between 
drains  operating  at  different  GRIs.  This  close  proximity  of  different  drains 
could  lead  to  a cross-rate  interference  problem  if  a LORAN-C  receiver  front 
end  were  not  properly  designed. 

It  is  additionally  possible  to  install  one's  own  LORAN  chain  by  tire  rental 
of  small  (1  kw)  portable  transmitters  similar  to  the  JOEY'-D  tactical  system. 

This  allows  a LORAN-C  drain  to  bo  set  up  anywhere  in  the  world  at  moderate  oost. 
Transmitters  are  available  from  Msgapulse,  Inc.  of  Bedford,  and  Irave  a range 
of  about  400  mi. 

Effect  of  Distance  on  SNR 

The  usefulness  of  a received  signal  is  expressed  as  the  signal-to- 
noise  ratio  (SN”R) . This  ratio  is  primarily  dr  pendent  on  distanoo  from  the 
transmitter , liowever  it  is  also  affected  by  other  factors  such  as: 


h. 
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Ficjurc  3 - I/OIVu-J-C  Positional  Repeatability  for  0.1  usee 

TD  Variation  on  East  Gaast  Grain 
(fron  Microlooic  Inc.  Onerators  f'.anual) 


— reflectivity  of  earth's  surface  (different  for  land, 
sea  water , fresh  water) 

— time  of  day  and  year 

— distance  from  the  equator 


Using  the  average  values  of  all  these  factors.  Figure  5 shews  SNR 
as  a function  of  distance. 


3.3.4  IORAN-C  System  Availability 

The  IORAN-C  system  has  maintain ad  a 993  availability  (Poland,  1973). 

If  scheduled  maintainance  is  not  included,  the  figure  beccr.cs  99.73.  Presently, 
nev  equipment  is  being  developed  which  will  permit  more  on-air  maintainance , 

which  will  improve  service  even  more. 


3.3.5  IORAN-C  Receiving  Equipment 

The  new  generation  of  IOFAN-C  receivers  currently  on  the  market  are 
ideally  suited  to  an  automatically  positioned  drifting  buoy  application.  Using 
microprocessors,  the  receivers  automatically  search  for  the  IORAN-C  signals, 
lock  on  to  the  tiiird  cycle  crossing,  and  compute  time  differences.  Such 
receivers  generally  monitor  the  SN*R  and  check  for  hi  ink  code  messages.  As 
described  earl  ier , LORAN-C  has  provisions  for  serti na  out  grovndwaves  from 
skywaves,  and  tills  is  done  by  all  fully  automatic  receivers.  Sena  manufacturers 
have  even  proviccd  a frequency  shift  keying  (ESP)  option  cn  the  output  such  that 
tiie  time  difference  words  acre  converted  to  signals  that  con  nxlulate  a higher 
frequency  data  return  link.  Tiiis  option  is  presently  available  on  both  Uie 
Teledyne  and  Intemav  micrcprocessor  LORAN-C  receivers. 


3. 3.5.1  Power  Requirements 

Power  consumption  for  automatic  receivers,  with  displays  eliminated, 
varies  fran  15  watts  to  38  watts.  In  either  case,  it  would  be  impractical  to 
leave  the  receiver  on  during  the  entire  experiment.  Eecause  a receiver  can 
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automatically  acquire  and  track  LORAN-C  signals,  it  is  possible  to  turn  on 
the  receiver  for  only  a short  time  for  each  sartple.  After  the  receiver  has 
searched,  settled,  and  computed  a pair  of  time  differences,  the  two  time 
difference  words  could  be  immediately  transmitted  or  stored  m a random 
access  memory  and  the  receiver  shut  down  until  the  next  sample  is  to  be  taken. 
Hie  total  power  consumption  for  only  the  acquisition  of  time  difference  infor- 
mation is  therefore  dependent  on  the  sampling  rate,  the  amount  of  power  drawn 
by  tiie  receiver,  as  well  as  the  amount  of  time  required  for  the  receiver  to 
determine  the  two  tine  difference  words.  This  is  called  search  and  settle 
time. 

Search  and  Settle  Time 

Search  and  settle  time  varies  with  -the  receiver  design  aid  is  de- 
pendent on  SNR.  Figure  6 shows  search  and  settle  time  as  a function  of  SNR 
for  four  leading  manufacturers  of  LORAN-C  receivers.  Hie  great  disparity  in 
receiver  performances,  seems  to  vary'  with  the  sampling  teciiniques  on  the 
front  end  and  efficiency  of  use  in  tne  microprocessor.  The  Nortlistar  G000  re- 
ceiver made  by  Digital  Marine  Electronics  in  Bedford,  Mass,  appears  to  be 
outstanding  on  this  point  plus  field  report's  on  this  svstem  are  glowing. 

The  Internav  IC-123  is  too  new  to  knew  much  about  it  and  thus  the  dashed  curve 
in  the  lew  SNR  region  of  Figure  6.  The  Telcdyne  701  data  derived  from  the 
manufacturer  are  felt  to  be  somewhat  optimistic  in  light  of  test  data  derived 
by  the  Coast  Guard  on  tills  unit.  Hie  Micrologic  ML-200  data  are  felt  to  bo 
realistic  upper  time  bounds  for  power  budget  estimates  for  the  case  of  a Gulf 
Stream  study  for  the  following  reasons: 

(1)  Tna  u.  S.  Coast  Guard  Electronics  Research  lab  in  Wildwood, 

New  Jersey,  whidi  has  been  testing  IORAN-C  receivers  which 
wore  marketed  up  to  about  1-year  ago,  indicate  that  reliable 
master  and  slave  station  lock-on  can  only  be  obtained  under 
law  SNR  conditions  by  the  longer  times  shown  in  Figure  6. 
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(2)  With  the  addition  of  2 new  East  Coast  LORAN-C  cliains 
in  1978,  the  cross-rate  interference  problem  between 
chains  of  different  GRIs  nay  rrake  search  plus  settle 
tines  greater  than  originally  measured  by  soire  receivers. 

The  cross-rate  interference  problem  may  not  be  present  in  tracking 
the  Kuroshio  (only  one  IDRAN-C  chain)  but  conservative  tine  estimates  must 
be  used  in  the  calculations. 

3. 3. 5. 2 Net  Power  Consumption 

The  amount  of  battery  energy  which  must  be  stored  to  pa  war  the 
receiver  over  the  entire  30 -day  experiment  nay  be  expressed  as: 

DS3PT 

E = 3600  (2) 

where 

E = energy  storage  required  (watt -hr) 

D = number  of  days  of  experiment  = 30  (nominal) 

S - no.  sanples/day 

3P  = power  drawn  by  receiver  (watts) 

Ts  = average  search  and  settle  time  (sec) 

It  can  be  seen  that  the  amount  of  batteries  required  will  be  directly  dependent 
on  receiver  power,  sampling  rate,  and  searcli  and  settle  time.  Thus,  both 
power  and  speed  of  search  and  settle  must  be  considered  in  order  to  determine 
how  much  energy  a receiver  will  require. 

Table  1 is  a sunnary  of  the  30-day  energy  requirements  as  well  as 
other  pertinent  parameters  for  four  leading  manufacturers  of  LORAN-C  receivers 
that  could  be  employed  in  a drifting  buoy  application.  The  details  of  the 
pev/er  budget  calculation,  and  the  sensitivity  to  sampling  interval,  are  given 
in  Appendix  A.  The  receiver  SNR  employed  in  the  calculations  (i.e.,  0 dB) 
is  that  of  a typical  night-time  SNR  (lower  than  daytime)  at  a distance  of 
approximately  1850  km  (1000  miles)  from  the  station.  Tliis  is  felt  to  be 
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\T-lonufacturer 

Digital  Marine 

and 

Model 

Micrologic 

Electronics 

Intcmav 

ML-200 

Northstar 

(SIMRAD) 

Teledyne 

Parameter  xvN^ 

6000 

i LC-123 

Model  701 

Power  Inunctions : 

Power  Drain 

15.25  w. 

33  w. 

'vlS  w. 

20.4  w. 

T @ 0 dB  SNR 
s 

^310  sec 

w60  sec. 

^240  sec. 

w90  sec. 

DC  powered 

yes 

yes 

yes 

yes 

Voltages 

12/2 4/3 2vdc 

10-40  vdc 

12/32  vdc 

12  vdc 

Front  End  Filters 

No. Fixed  Notches 

4 

0 

None 

0 

Variable  Notches 

4 

4 

None 

2 

Rockwell 

Microprocessor 

Intel  4004 

Intel  3030 

Intel  8030 

pps-4 

Logic 

TTL 

TTL 

La//  pwr 

TIL 

Schottky 

Slaves  Tracked 

5 

5 

4 

4 

Blink  Code 

yes 

yes 

no 

no 

SNR  Tracking 

yes 

yes 

no 

yes 

Memory  Avail. 

extra  cost 

extra  cost 

extra  cost 

yes  (option.) 

FSK  Output 

no 

no 

yes  (option) 

yes  (option . ) 

Fisiiing 

Ccmrercial 

New  system 

Considerable 

Field  Experience 

boats 

boats 

very  good 
results 

very  good 

excellent  | 

little  1 

experience 

Price:  (*) 

Single: 

"v.$2600.  (GSA) 

$3,200  (GSA) 

w$2450. (GSA) 

w$2900. (GSA) 

Lot  qty: (10+) 

^$1925. (GSA) 

$3200. (GSA) 

$2340. (GSA) 

n$2900. (GSA) 

30-Day  Energy 
Budget  (1  pt/hr) 

946. 

456. 

864. 

367. 

(watt-hrs) 


* Includes  antenna  coupler  ("-$200. ) 

Table  1 - Surmviry  Comparison  of  leading  LOPAN-C 

Receivers  Adaptable  to  Drifting  Buoy  Applications 
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representative  of  a maximum  distance  found  while  tracking  a western  boundary 
current  such  as  the  Gulf  Stream.  The  large  deviations  in  the  estimated 
search  (i.e.f  signal  acquisition)  and  settle  (i.e.,  cycle  selection)  times 
of  the  various  receivers  seems  to  arise  from  design  differences  in  the  re- 
ceiver front  end.  It  appears  tint  both  the  horthstar  6000  and  the  Intemav 
IC-123  designs  have  placed  a special  emphasis  on  tliis  aspect  by  achieving 
high  signal  power  gains  while  maintaining  linearity.  In  so  doing  it  appears 
that  more  power  may  be  dissipated  in  this  portion  of  their  units  than  with 
other  designs.  This  front  end  signal  conditioning  may  be  especially  important 
on  a drifting  buoy  that  would  be  retransmitting  an  HF  or  Ui.IF  signal  through 
aii  adjacent  antenna  tint  may  create  a "noisy"  environment. 

The  estimated  prices  shewn  in  Table  1 were  derived  in  a somewhat 
different  manner  for  the  Micrologic  ML-200  system  than  for  the-  other  units. 
The  ML-200  list  price  is  based  on  a breakdown  of  prices  for  the  key  system 
elements:  digital  computer  board  ($1566.),  LORAlv-C  receiver  ($330.),  antenna 
coupler  ($192.),  and  antenna  (v$100.).  The  prices  of  the  other  three  units 
are  estimates  for  a full  system,  including  case,  display  lights,  buttons,  and 
antenna  coupler,  because  that  is  the  simplest  manner  by  wliich  tiiey  can  sell 

systems  and  estimate  price.  The  equivalent  G3A  list  price  on  the  Micro- 
logic ML-200  is  $3775.  ($3104  for  10  or  more).  It  would  be  a task  for  the 
system  fabricator  to  remove  the  case  and  disconnect  the  lights.  Perhaps 
this  pricing  aspect  would  change  favorably  for  such  systems  as  the  Northstar 
6000  if  a more  serious  price  quotation  were  sought  on  a quantity  purcliase. 

It  should  also  bo  mentioned  that  approximately  $100.  should  be  added  to  the 
costs  given  in  Table  1 to  account  for  an  S'  LCPAh-C  whip  antenna. 

Many  other  manufacturers'  catalogs  were  consulted  but,  at  the  time  of 
this  writing,  their  products  were  deemed  to  be  not  fully  automatic  enough  for 
buoy  operation.  For  example,  some  units  required  visual  LORAN-C  cycle  notch- 
ing with  a cathode  ray  tube;  while  others  did  not  contain  a microprocessor  for 
acquiring  and  tracking  the  LORAN-C  signal.  In  the  future,  if  the  design 
suggested  were  to  be  inplemented,  the  available  units  should  be  resurveyed 
because  of  the  extremely  competitive  and  rapidly-clianging  market  for  such 
systems.  A few  factors  make  for  both  a technical  and  marketing  volatility  for 
LORAN-C  systems: 
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(1)  Hie  advent  of  microprocessors  and  their  lowering 
price  and  expanding  capabilities. 

(2)  Hie  expanded  capability  on  the  LOPAN-C  system  will  bring 
in  more  land-based  users. 

(3)  The  impending  shut-down  of  the  LORAh-A  navigation  system 
in  approximately  1981  forcing  fishermen  and  small  teat 
owners  to  use  LOPAN-C. 

(4)  The  new  200-mile  fishing  limit  and  its  patrol  - possibly 
necessitating  foreign  vessels  to  employ  automatic 
LOPAN-C  to  ba  assured  of  compliance. 


It  does  not  appear  as  though  any  LORAN-C  receiver  manufacturer  will 
build  a new  receiver  employing  low  power  CDS  (complimentary  metal  oxide 
semiconductor)  logic.  There  are  three  very  goad  reasons  for  not  doing  such  : 


(1) 


CDS  logic,  including  even  the  more  reliable  ceramic 
encapsulated  version,  has  not  proven  to  bo  sufficiently 


relic 
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average  LORAN-C  receiver  would  be  subjected  at  sea.  The 
average  market,  which  is  being  addressed  by  the  manufacturers, 
is  that  of  fishing  vessels  and  larger  siiips  on  which  the 
receiver  is  installed  in  the  pilot  house  - subjected  to  the 
hazardous  environment  of  humidity,  vibration,  and  salt  air. 


(2)  Tnere  is  no  strong,  compelling  reason  to  reduce  the  power 
consumption  of  the  receiver  from  a range  of  40  to  75  watts 
(including  display  lights)  to  approximately  20  or  30  watts 
(again  including  lights)  when  vessel  power  in  these  ranges 
is  readily  available. 

(3)  There  are  not  as  many  desirable  design  components  available 
employing  CDS  logic. 


If  a CMOS  version  of  the  stripped-down  receivers  described  in  Table  1 
(i.e.,  no  display  lights)  were  implemented,  the  power  dissipations  would  drop 
to  the  3-5  watt  range  - consum'd  primarily  in  the  front  and  receiver  and 
antenna  coupler. 
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3.4 


Navy  Navigation  Satellite  System  (NNSS)  or  TRANSIT 


The  Navy  Navigation  Satellite  System  (NNSS) , often  referred  to  as  the 
TRANSIT  System,  lias  been  in  existence  since  1964 . It  consists  of  5 polar- 
orbitting  satellites  (1.8  hour  orbittal  period)  which  permit  a position  fix 
approximately  every'  hour  and  one-iialf  to  accuracy  of  .04  km  to  2 Ion  de- 
pending on  the  receiver  sophistication  and  knowledge  of  receiver  velocity. 

In  order  to  obtain  the  highest  accuracy  a user,  such  as  a large  ocean-going 
ship,  generally  uses  a receiver  which  costs  over  $20,000.  It  has  been  sug- 
gested by  Vfcsterfield  (1972)  tliat  a simplified  TRANSIT  receiver  be  installed 
aboard  a drifting  buoy  for  the  acquisition  of  trajectory'  data  only,  anywhere 
in  the  v.orld.  Such  a receiver  would  cost  approximately  $5,000  (in  quantity), 
give  an  accuracy  of  the  order  of  2 km,  a fix  approximately  every  1.5  hours, 
and  require  approximately  660  watt-hours  of  electrical  energy  over  a 30  day 
mission.  Such  power  consunption  would  necessitate  approximately  40  pounds  of 
alkaline  cells.  The  suggested  system  provides  no  capability  for  sending 
additional  ancillary  data  such  as  a drogue  indicator,  wind  spaed,  or  barometric 
pressure.  Because  many’  of  the  above  factors  do  not  neat  the  system  objectives 
the  TRANSIT  system  was  rejected  at  present.  It  is,  hewever,  possible  tlaat  in 
the  future  a low  power-consuming  receiver  ray'  be  efficiently  manufactured  in 
quantity  with  adequate  accuracy  such  that  it  may  become  more  competitive  with 
other  systems  described.  In  light  of  the  fact  that  the  Global  Positioning 
System  (GTS)  will  replace  the  function  of  TRANSIT  in  the  1930's  it  is  felt  that 
such  a TRANSIT  receiver  will  not  be  made. 
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3.5  NAVSTAR  Global  Positioning  System  (GPS) 

The  Global  Positioning  System  (GPS)  is  a high  accuracy  satellite 
positioning  system  being  designed  primarily  for  the  U.  S.  Military  under 
Air  Force  sponsorship.  Until  recently  both  the  Navy  and  Air  Force  were 
pursuing  designs  of  similar  systems  - the  Navy  system  being  called  Tarnation 
was  destined  to  be  the  successor  to  TRANSIT.  For  efficiency  all  efforts 
have  recently  been  channeled  into  the  GPS  effort  only.  The  GPS  system  is 
expected  to  exhibit  an  overall  positional  accuracy  (including  all  error 
sources)  of  appronirra t e ly  10-15  meters  independent  of  velxicle  velocity  with 
positions  provided  on  command.  TRANSIT , cn  the  other  hand,  is  inherently 
velocity  sensitive  (unless  very  costly  equipment  is  used)  and  only  provides 
fines  on  a periodic  basis  wiiile  requiring  a total  20-minute  pass  in  order  to 
acquire  a fix. 

The  GPS  system  will  eventually  consist  of  20  polor-orbitting  satellites 
(12  hour  period)  such  that  -}  satellites  are  always  in  view  at  any  time.  The 
present  schedule  colls  for  a 2-dimcnsional  GPS  capability  to  be  ready  in  1980 
(good  for  surface  oceanography)  and  a 3 -d in>_ 1 i s ions  L system  ready  by  1381 
(Dennis,  1373)  • The  satellites  wij.x  Li  e euxici  - tequeiiCies  — one  Sig- 

nal coded  in  a jam-resistant  manner  for  the  military’  oxi  another  coded  in  a 
well-defined  manner  for  all  types  of  use.  It  appears  tint  tiac  accuracy'  of  the 
non-military  channel  may  be  poorer  than  that  described  already  but  still  more 
than  adequate  for  the  needs  of  drifting  buoys . 

At  present  it  appears  that  one  manufacturer  ray  try  to  come  out  with  a 
"Spartan"  GPS  receiver  which  they  will  attempt  to  market  for  approximately 
$2000.  There  is  no  word  yet  on  power  drained  from  such  a unit.  lastly, 
although  the  GPS  system  may'  be  available  for  general  usage  there  ray  be  a use 
cliarge  associated  with  it  wiiich  should  be  factored  into  system  cost  analysis. 

In  summary  it  appears  that  the  GPS  system  has  many  features  which  are 
appealing  for  a drifting  buoy  positioning  system  of  the  future.  Because  of 
its  schedule  and  unknowns  associated  with  buoy  system  costs  and  pcwer,  the 
system  cannot  be  fully  examined  at  present. 
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3.6  Other  Positioning  Systems 


Other  drifting  buoy  remote  positioning  systems  are  available  but  are 
generally  of  such  low  accuracy  or  limited  range  that  they  cure  inapplicable 
for  tire  study  of  higher  frequency  motions  of  western  boundary'  currents.  For 
example,  Kirwan  and  McNally  (1975)  tracked  drifting  buoys  in  tire  North  Pacific 
with  over-the-horizon  (OfH)  radar  receiving  signals  from  a transmitter  on 
the  buoy.  The  buoys  were  tracked  to  ranges  of  approximately  3000  km.  with  a 
positional  accuracy  of  approximately  -8  km.  In  addition,  Whelan  et  al  (1975) 
of  tire  ITT  Electro  Physics  lab  designed,  built,  and  marketed  a law  cost 
drifting  buoy  positioning  system  using  an  I IF  transmitter.  A - 1/2°  azimuth 
resolution  was  claimed  on  the  12?  direction-finding  system.  Tire  limited  range 
and  accuracy'  of  such  a system  (a 300  Ion,  max.),  problems  with  night-time  IIF 
radio  interference,  plus  tire  cost  of  manning  tv.o  shore  stations  for  signal 
triangulation,  were  factors  making  tire  system  unuseful  for  tire  purposes  des- 
cribed. 
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4 . 0  TELET’TTRY  OF  BUOY  POSITION  AND  SENSOR  DATA 


The  telemetry  of  data  frcm  the  buoy  is  the  major  problem  associated 
with  the  given  design.  There  are  only  two  means  for  the  transmittal  of  buoy 
data  - high  frequency  (IIF)  radio  signals  employing  frequency  shift  keying 
(FSK)  or  the  use  of  satellite  transmitters.  Hiare  is  a choice  of  two  rela- 
tively law  cost  satellite  systems  for  the  expected  buoy  deployments  of  1979- 
1930,  the  TIROS-N  or  the  GOES  satellite.  Each  system  of  data  return  has  its 
positive  and  negative  aspects  and  will  be  described  in  this  section. 

4.1  Satellite  Tblenctry 

The  TIROS-N  and  GOES  satellites  are  the  cnly  real  candidate  satellites  for 
drifting  buoy  applications  because  of  their  low  system  cost  and  power  consump- 
tion. Other  systems,  available  or  projected  for  the  transmittal  of  data,  are 
possible  candidates,  but  tiieir  system  or  usage  costs  (or  both)  are  too  high. 
These  systems  will  be  mentioned  in  Section  4.1.3. 

4.1.1  TIROS-N  Satellite  Telemetry 

The  TIROS-N  Buoy  Transmit  Terminal  (BIT)  is  a low  cost  transmitter 
(v$16C0.  with  antenna,  single  qty)  that  employs  a reliable  telemetry  system 
with  considerable  experience  derived  from,  its  progenitors,  the  ETi '/NIMBUS- 6 
system  and  the  IRLS  (Jhterrogauion,  Recording  and  location  Subsystem)  aboard 
tiie  NIMBUS  II  and  III  satellites.  The  TIROS-N  satellite  should  be  launched 
and  available  for  use  in  mid-1978.  NASA,  Goddard  Space  Flight  Center,  who 
operates  the  NIMBUS  and  TIROS  satellites/  premises  cnly  one  position  aid 
telemetry  link  per  day,  from  each  satellite  at  * 5km  (i.e. , 3 a data)  accuracy, 
but  investigators  have  been  deriving  data  at  frequencies  of  2 to  8 times  per 
day  from  a single  NIM3US-6  satellite.  These  positions  are  accurate  to  approx- 
imately * 1-2  km  with  repeatabilities  of  less  than  - 1 km  (i.e.,  3 o limits). 
There  will  be  2 sun-synciironous , polar- orbit ting  TIROS-N  satellites  with  a 
108-minute  orbital  period.  Etch  satellite  should  be  visible  at  least  twice 
per  day.  Because  only  one  satellite  will  be  operable  for  a while,  awaiting 
the  launch  of  the  second,  and  because  one  satellite  may  become  inoperable  at 
a future  date,  the  TIROS-N  data  telemetry  link  will  be  conservatively  based 
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on  a single  satellite.  If  the  2 satellite  systems  were  analyzed  the  design 
would  vary  little. 

It  will  be  assumed  that  a half-redundant  data  point,  with  the  word 
structure  for  each  data  point  as  shewn  in  Appendix  A,  will  be  radioed  up  to 
the  satellite  every  12  hours.  The  4 least  significant  digits  of  a LORAN-C 
ID  (including  tenths  of  microseconds)  are  coded  as  12-bit  words. 

A strawraan  LORAN-C/TIROS-N  hybrid  buoy  positioning  system  is  describe 
in  Figure  7.  A typical  power  budget  and  cost  for  such  a system  is  given  in 
Table  2 based  on  the  Micrologic  ML-200  microprocessor  LORAN-C  receiver  described 
in  Table  1.  This  receiver  is  chosen  because  of  its  adequate  performance  and 
low  cost.  The  key  features  of  the  system  which  make  it  more  complex  and  yet 
very  flexible  are  the  system  sequencing  logic  and  on-board  memory  - both  of 
which  are  necessarily  powered  all  of  the  time.  Because  of  their  lew  power  drains, 
the  temperature  and  bar  cm?  uric  pressure  sensors  as  well  as  the  drogue  indicator 
are  assumed  to  be  powered  up  all  of  the  time.  The  memory  is  assured  to  be 
necessary  because  the  LORAN-C  ID's  must  be  stored  between  periodic  position 
fixes,  awaiting  telemetry  to  the  satellite.  The  high  power  requirements  of  the 
present  microprocessor  LORAN-C  receivers  require  that,  between  position  fixes, 
the  unit  be  powered  down.  If  such  were  not  the  case  the  receiver  energy  consump- 
tion would  increase  15-fold  from  its  already  burdensome  946  watt-hours  for 
hourly  fixes. 

Wherever  possible  G-DS  logic  is  employed  in  order  to  conserve  power. 
Special  care  should  be  taken  in  order  to  employ  high  reliability,  ceramic- 
encapsulatod  components.  The  estimated  prices  reflect  single  unit  purchases 
(except  for  the  LORAN-C  receiver)  after  developmental  or  design  work  is  com- 
pleted on  the  timing  and  memory  circuitry.  Gare  will  be  taken  not  to  require 
alteration  of  the  manufacturer's  microprocessor  functions  and  programs.  Certain 
elements  in  Figure  7 are  eliminated  for  simplicity.  For  example,  it  is  assumed 
that  tiis  signal  frem  the  existing  of f-tho-shclf  microprocessor  receiver  is 
readily  inputted  to  a random  access  memory  (RAM)  or  to  a simple  first- in,  first- 
out  buffer  memory.  These  signals  will  in  most  cases  have  to  first  be  converted 
to  the  proper  form  (parailcl-to-scrial  conversion,  etc.)  for  acceptance  by  the 
memory.  In  like  manner,  the  signal  from  the  memory  to  the  satellite  transmitter 
nay  require  changes  in  order  to  interface  properly.  These  additions  add  very 
little  to  the  system  cost,  because  components  are  inexpensive,  but  do  complicate 
the  engineering  somewhat . The  next  section  well  indicate  a design  that  is  much 
simpler  but  more  expensive. 
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Component 

Current 

Drain 

(ma) 

Voltage 

(Volts) 

Power 

Drain 

(Watts) 

Energy/30 
days  @ 

1 fix/hr 
(v.’att-hrs) 

Component 

price 

after 

Development 

DORAN- C 
y-proc . revr 
and  antenna 

1270 

Assume 
+ 12 

15.25 

946 

v$2100.- 

Clock*  & 

Sequencing 

Functions 

10  ma. 

+5 

.05 

36. 

1 

o 

CD 

</> 

CMOS* 

Memory 

.000020 

+5 

.0001 

.07 

$10.- 

Water 

Temperature 

.00024 

+12 

.003 

2.2 

^$100.- 

^Jxxro:TC  tr  i c+ 

1.0 

+12 

.012 

8.6 

O/$500.~ 

Drogue 

Indicator 

nil 

+12 

nil 

nil 

v$100.- 

TIROS-N 

mm  + 

Antenna 

550 

+12 

6.6  (max) 

83.9 

$1600.- 

TOTALS 

1831 

- 

21.92 

f ; 

1076. S 

K 

$4490.- 

Note:-  Estimated  power  drains  include  pov.cr  supply  voltage  regulation. 
-*Asterisks  indicate  required  design  work  - others  off-the-shelf. 
-+Dunkcr-Rtrno  (NCAR)  barometric  pressure  all  assumed. 

4For  Itorthstar  C000  unit:  586  v.— Hrs,  $5690.  (tot.) 


Table  2 - Typical  Power  Budget  and  Unit  Cost  for  a 

IORAN-C/TIPOS-N  Hybrid  Drifting  Buoy 
Positioning  System 


24 


4.1.2  GOES  Satellite  Telenotry 


r 


Because  the  GOES  satellites  are  synchronous  over  the  equator  (75°  w. 
long,  100°  w.  long,  and  135°  w.  long)  a drifting  buoy  positioning  itself  by 
LORAN-C  can  immediately  get  access  to  GOES  for  data  transmittal.  In  this 
way  the  system  does  not  require  an  on-board  memory  capability.  The  design 
problem  is  reduced  to  mainly  that  of  a clock/sequcncing  function  for  the  LORAN- 
C microprocessor  and  proper  data  handling  (see  Figure  7).  The  power  budget 
and  costs  for  such  a system  (single  qty)  are  summarized  in  Table  3.  This 
nominal  power  budget  assumes  a pair  of  12-bit  JjORAN-C  TD's;  one  12-bit 
measurement  of  barometric  pressure  and  water  temperature,  and  a 4-bit  drogue 
indicator  word  - totalling  52  bits  of  data  transmitted  every  hour.  The  timed 
acquisition  of  the  'ID's  is  governed  by  the  clock  and  sequencing  circuitry. 

As  soon  as  those  data  are  ready  (after  search  and  settle  time)  they  are  fed 
directly  to  the  GOES  in  a BCD  format.  They  arc  immediately  transmitted  along 
with  the  P Dtv-q,  T , , and  drogue  indicator;  which  are  instantaneously 

sampled  and  digitized  within  the  GOES  transmitter. 

The  GOES  transmitter  telemeters  up  to  16  words  per  transmission  at 
ICO  bits/scc  in  the  following  general  format ■ 

GOES  Timing  Format-  (hourly  Data  Transmissions) 

Carrier  Enable: 

Minch-Coded  bit  sync 
15-bit  (mux. ) sequence 
31-bit  cor.Tiund  word 
52-bit  data  word 

Tetal  Time:  8.48  sec 

The  transmitter  consumes  .05  watts  (quiescent)  between  transmissions  and  40- 
watts  during  the  above  estimated  transmission  period.  This  leads  to  a 30-day 
energy  budget  of  108.8  watt-hours  as  shown  in  Table  3.  The  resulting  total 
energy  requirement  for  the  IORAE-C/GOES  hybrid  system  is  estimated  to  be 
857.7  watt-hours  - requiring  approximately  45-pounds  of  alkaline  cells. 


5 . 0 sec 
2 . 5 see 
.15  sec 
.31  see 
.52  sec 
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Component 

Current 

Drain 

(rra) 

Voltage 
(Volts) , 
DC 

Power 

Drain 

(Watts) 

Energy/' 30 
days  @ 

1 fix/)  IT 
(watt-hrs ) 

Component 

Price 

after 

Development 

LORAN-C 
-proc.  revr 
and  antenna 

1270 

12 

15.25 

702 

%$2100.- 

Clock  and 
Sequencing 
Functions 

10  ma 

+5 

.05 

36.0 

$80. 

Q-DS 

Memory 

.00002 

+5 

.0001 

.07 

$10 

Water 

Tenperature 

.00024 

+12 

.003 

2.2 

a/$100 

^barometric 

1.0 

+12 

.012 

8.6 

^$500. 

Drogue 

Indicator 

mil 

+12 

mil 

nil 

a.$100. 

GOeS  ' 
Satellite 
Transmitter 
+ Antenna 

3300  (max) 

.04 

(quiescent) 

+12 

40  (max) 

.05 

(quiescent ) 

103.8 

$3300.- 

TOTALS: 

4531  (max) 

- 

55.3  (max) 

857.7 

$6190.- 

+ Assume  Bunkcr-Ramo  (XCAR)  p,  sensor 

1 boro 

£ Assuming  52  bits/!ir  (i.e.,  2 TD's,  p,  T,  drogue  indicator) 


Table  3 - Typical  Power  Budaet  and  Unit  Cost 

for  a I ORAX-C/GOF.S  Hybrid  Drifting  Buoy 
Positioning  System 
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The  GOES  energy  estimate  is  very  co-par  able  to  that  of  the  LORAN-C/ 
TIRDS-N  hybrid  system  wiiich  might  require  1C7G.8  watt-hours . The  GOES  system 
would,  however,  cost  approximately  $1700.  more  per  buoy  than  the  TIROS-N 
hybrid,  but  engineering  development  costs  associated  with  implementing  the 
TIROS-N  iTorory  would  somcwiiat  offset  this  cost  difference.  The  original 
proposal  called  for  the  deployment  of  six  expendable  drifting  buoys  every 
three  months  for  two  years  or  a total  of  43  buoys . At  a differential  cost  per 
buoy  of  $1700.  this  amounts  to  a total  cost  difference  of  $31,600.  over  the 
2 -year  life  of  the  experiment . If  the  IORAN-C/TIROS-N  system  shewn  in  Figure 
7,  were  designed  and  implemented  ‘die  memory  costs  shown  in  Table  2 are  negli- 
gible (i.e.,  u$10.  per  buoy).  The  engineering  and  test  costs  associated  with 
implementing  the  memory  and  timing  interfaces  to  the  TIROS-N  transmitter  are 
conservatively  estimated  to  bo  less  than  $20,000.  The  overall  battery-pack 
costs  and  reliability  of  each  system  in  terms  of  longevity  and  data  quality 
are  felt  to  be  equivalent.  Therefore,  it  is  estimated  that  the  IORAN-C/TIROS-N 
hybrid  design  would  cost  approximately  $60,000.  less  than  the  IORAN-C/ GOES 
system  over  the  2-yoar  life  of  the  planned  experiment. 


4.1.3  Other  Satellite  Data  or  Telemetry  links 

There  are  numerous  other  satellites  which  are  useful  for  either  military 
purposes,  earth  monitoring  proaroms,  or  data  and  voice  corrrranicaticns  from  sea. 
This  section  will  briefly  describe  seme  of  the  most  frequently-mentioned  systems 
to  varying  levels  of  detail. 

MARI SAT 

A communications  satellite  operated  by  Comsat  General  Corporation  for 
the  U.  S.  Navy  and  private  use.  It  provides  ship-to-shore  and  shore-to-ship 
telegraph  and  telephony  service  working  in  the  GHz  (private  use)  aid  UliF 
(U.  S.  Navy)  frequencies,  It  requires  a 1 . 2-meter  diameter  gimballed  (sta- 
bilized) antenna.  Systems  generally  find  usage  on  offshore  oil  rigs  and  explor 
ation  vessels.  Costs  for  purchase,  rental,  and  use  are  high  compared  to  the 
guidelines  in  this  report. 


J 
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GEOS 


A £Oodetic  earth  observation  satellite  with  no  communication  or 
location  capabilities.  This  typo  of  satellite  is  still  being  launched  by 
the  multi-nation  European  Space  Agency. 

LANOSAT  ( formerly  ECTS) 

LAXDSAT  is  tiie  new  name  for  the  NASA  Earth  Resource  Technology  Satellites 
(EKTS)  which  contain  multi-spectral  scanners  for  earth  monitoring.  There  are 
no  on-board  communication  or  location  capabilities  on  LANDSAT. 


SBASAT-A 

SEASAT-A  is  part  of  tlio  NASA  applications  satellite  program.  Ib  be 
launched  in  mid  1973,  it  will  carry  a radar  altimeter,  capable  of  measuring  sea 
surface  topography  (including  currents  and  tides);  a radar  scatterometer  in 
order  to  measure  surface  winds  by  their  effect  on  capillary  waves;  a scanning 
multi-frequency  microwave  radiometer  for  the  measurement  of  sea  surface  tem- 
perature; and  a synthetic  aperture  radar,  capable  of  providing  data  on  ocean 
waves,  coastal  regions,  and  sea  ice  (personal  communication) . 

Estimates  of  the  locations  of  currents  can  be  inferred  through  measure- 
ments of  sea  surface  topography  and  temperature.  The  temperature  resolution  of 
the  multi-frequency  radiometer  will  vary'  from  0.34  to  1.09°C,  depending  on  the 
frequency  used.  This  temperature  resolution  represents  an  average  over  a foot- 
print whose  smallest  size  is  21  x 14  km.  The  temperature  resolution  veuld 
be  potentially  adequate  for  the  study  of  boundary  currents  but  the  spatial  reso- 
lution is  too  coarse.  An  additional  modest  visual  and  infrared  imaging  instru- 
ment will  be  aboard  SEASAT-A  but  it  will  only  be  capable  of  a 1°C  resolution 
with  a surface  resolution  of  9 km  - again  too  coarse. 

The  performance  and  data  from  SEASAT-A  should  be  closely  followed  as  a 
strong  adjunct  to  neor-surfacc  studies  of  not  only  western  boundary  currents 
but  such  areas  as  rings  and  internal  waves . The  capability  to  remotely  monitor 
tenperature,  sea  state,  and  wind  field  data  may  provide  the  necessary  data  base 
for  more  accurately  assessing  the  quality  of  Lagrangian  drifter  data  (i.e., 
slippage  estimation) . 
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ATS-5 


The  NASA  Applications  Technology  Satellites  allow  the  transmittal 
of  voice,  teletype,  facsimile,  arid  coded  data  to  shore  via  geostationary 
satellites.  Positions  to  an  average  accuracy  of  approximately  2.0  km  are 
also  available.  The  systems  are  finding  increased  use  aboard  ocean  research 
vessels  to  effectively  link  the  vessel  to  shore-based  computers.  The 
system  v.ould  be  inappropriate  for  a drifting  buoy  application  because  it 
requires  a directional  or  steered  antenna  on  the  buoy. 
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4.2  High  Frequency  (HF)  Data  Telemetry 

The  feasibility  of  using  frequency  shift  keying  (FSK)  single  sideband 
(SSB)  high  frequency  telemetry'  Iras  been  examined.  High  frequency  telemetry 
enables  the  experimenter  to  know  tire  buoy's  position  immediately  after  the 
I0RAN-C  cordinates  have  been  transmitted,  without  tire  delay  and  possible  cost 
required  for  satellite  services  or  the  costs  and  problems  involved  in  using  a 
privately-purchased  TIROS-N  ground  station  interrogation  unit.  It  is  not 
possible  to  use  tire  LORAN-C  pulses  to  directly  modulate  an  HF  carrier  because 
of  HF  bandwidth  limitations.  It  might  be  possible  to  use  a hard-limited 
LORAL’ -C  signal  (i.e.,  square-wave)  for  modulation,  but  timing  delays  in  IIF 
propagation  would  seriously  degrade  the  IORAL-C  accuracy. 

A buoy  mounted  HF  transmitter  is  cheaper  and  mere  flexible  than  a satellite 
system,  but  monitoring  stations  would  have  to  be  maintained  and  reliability 
would  not  be  as  good  as  with  a satellite  system  (Livingston  et  al,  1977). 

The  system  considered  here  is  very  simple  in  order  to  remain  within  the 
economic  constraints.  A m.d  ti- freon an  cv  system  such  as  that  described  by 
Livingston  et  al  (1977)  or  that  in  the  report  by  the  Intergovernmental  Oceano- 
graphic Commission  (1967)  could  bo  employed  in  order  to  increase  reliability, 
but  eysterc  cocts  und  pewex  ccnsixi^tion  v.oulcl  incro^.so  do'^irm  nuiclc^— 

lines.  A single  frequency  system,  besides  being  simple,  avoids  problems  of 
antenna  detuning  and  radiating  inefficiency,  while  employing  an  on-board  timer, 
data  memory  and  ionospheric  predictions  in  order  to  optimize  the  quality  of 
transmission . 

In  order  to  increase  the  data  reliability  of  a single  frequency  system 
it  will  be  assumed  that  transmissions  are  made  at  local  noon  and  midnight.  Each 
transmission  will  contain  the  data  from  the  previous  24  hours,  with  each  trans- 
mission repeated  three  times.  This  format  results  in  each  data  point  being 
transmitted  six  times  - three  times  at  noon  and  three  tirr.es  at  midnight.  Later 
discussion  will  show  that  around  midnight  a particular  frequency  choice  will 
lead  to  the  maximum  probability  of  reliable  skvwave  data  transmittal  on  a year- 
round  basis  while  maximizing  noon- time  ground-wave  propagation.  This  frequency 
could  bo  changed  to  optimize  predictions  for  a particular  time  of  year  - but 
only  within  allotted  frequency  bands.  The  transmissions  from  the  system  should 
be  self-timed  with  an  output  power  of  100  watts  to  an  8-foot  whip  antenna.  The 
nominal  data  rate  is  conservatively  assumed  to  be  100  bits/sec. 
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Supplying  power  to  the  transmitter  is  a primary  consideration.  Altliough 
the  transmitter  draws  a large  amount  of  power,  transmissions  are  short  enough 
such  that  the  total  power  consumed  over  a 30-day  experiment  is  not  impractical. 
'Hie  total  amount  of  time  that  the  transmitter  must  be  on  is  determined  by  the 
total  amount  of  data  to  be  transmitted  each  12  hours.  Based  on  a nominal  one 
position  sanplc/hour,  (i.e.,  2,  12-bit  TD's)  with  water  temperature,  barometric 
pressure,  and  the  drogue  indicator  being  sampled  every  12  hours,  the  data  format 
shewn  in  Table  4 is  assumed  (see  also  Appendix  B) . 

Table  4 

HF  Transmitter  Data  Word  Structure 


Parameter 


Bits/V’ord 


Max.  Words,/12  hr 


Drogue  Indicator 


data  bits  rccordcd/12  hrs.  : 
data  bits  transmit tcd/12  hrs 


Total  Bits 


LORAN-C,  TD1 

12 

12 

144 

LORAM-C,  TD2 

12 

12 

144 

Water  Temperature 

12 

1 

12 

^baro 

12 

1 

12 

316 

632 


Assuming  that  632  data  bits  must  be  transmitted  every  12-hours  results  in  an 
HF  transmission  time  computed  as  follows: 


Transmit  Format 


(1) 

Warm-up 

5.0 

sec 

(2) 

Carrier  enable  and  bit 
sync  (from  Picquenard  (1974) 

7.5 

sec 

(3) 

632  bits/lOO  bits/sec 

6.3 

see 

Total  Transmitter  On-Time: 

18.8 

sec 
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It  is  assumed  that  six  such  transmissions  are  made  each  day  with  approximately 
100-watts  of  power  to  the  antenna  while  consuming  approximately  200  watts 
total.  As  a power  safety  margin  it  is  assumed  that  during  the  warm-up  time 
the  transmitter  is  also  consuming  200  watts.  During  the  "off"  time  there  is 
no  standby  power  consumption  because  the  system  is  assumed  to  be  fully  acti- 
vated by  a low  power  clock  in  the  same  manner  as  the  LOPAN-C  receiver . The 
above  power  profile  leads  to  a total  30-day  energy  consumption  computed  as 
follows : 

lEjjp  = 30  days  (6  transnyday)  (18.8  sec/trans) 

or: 

E = 188  watt-liours 

Ihis  amount  of  energy  could  be  provided  by  approximately  10-pounds  of  alkaline 
batteries. 


Frequency 

A frequency  in  the  4.1  mliz  oceanographic  band  (i.c.,  4.1625  to  4.1560  mHz) 
was  chosen  for  this  analysis.  This  choice  was  made  to  allow  a large  portion  of 
the  Gulf  Stream  to  be  overed  using  reliable  ground-wave  propagation.  While 
this  choice  of  frequency  eliminates  the  possibility  of  long  range  skywave 
propagation  during  the  day,  adequate  skywave  propagation  can  be  expected  at  night 


Signal  Strength 

The  foremost  problem  associated  with  high  frequency  telemetry  is  pre- 
dicting the  signal-to-noise  ratio  for  different  times  of  day,  year,  and  sunspot 
cycle.  High  frequency  waves  are  propagated  by  many  modes,  and  the  background 
noise  is  also  liighLy  variable.  If  frequency  shift  keying  (FSK)  modulation  is 
used,  a shore-based  rrcoivxl  signal-to-noisc  ratio  of  10  db  would  be  required 
according  to  Picquenard  (1974). 

4.2.1  Ground  Wave  IE  Telemetry  Link 

The  most  reliable  propagation  mode  available  is  the  ground  wave.  No 
ionospheric  reflection  is  required,  so  this  mode  is  less  subject  to  the  cyclic 
variations  which  dominate  the  other  modes. 
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The  most  likely  frequency  for  ground-wave  use  is  4 mHz.  The  ITT  hand- 
book entitled  "Reference  Data  for  Radio  Engineers"  indicates  that  atmospheric 
noise  increases  by  about  55  db/decade  below  4 mHz.  Above  4 mHz,  ground-wave 
propagation  losses  increase  by  about  23  db/decade.  Additionally,  above  5 mHz 
scattering  and  ground-wave/ skywave  interaction  degrade  ground-wave  propagation. 

The  ground-wave  field  over  sea  water,  for  1 lev’  of  radiated  power, 
appears  in  Figure  3,  taken  from  Picqucnard  (1974).  The  atmospheric  noise  level 
was  also  conputed  from  figures  in  Picqucnard  's  book.  Figure  9 is  an  example  of 
such  a table.  The  noise  data,  nominally  presented  in  db  above  1 watt,  have  been 
converted  to  db  above  lpv/m  (see  Picqucnard,  p.212)  for  compatibility  with  the 
ground-wave  field  data.  Fluctuations  and  uncertainty  margins  are  included  in 
order  to  make  the  noise  data  valid  for  951  of  the  time  at  the  worst  time  of  year. 
Calculations  are  shown  in  Appendix  B. 

Antenna  gain  is  the  remaining  factor  needed  in  order  to  predict  receiver 
signal-bo-noise-ratios  for  the  ground-wave  mode.  The  following  conservative 
assumptions  will  be  made: 

Receiving  antenna  gain  0 db 
Transmitting  " " -15  db 

Additionally,  10  log  ~ ”10  db  must  be  added  to  the  curves  for 

field  strength,  as  those  curves  are  for  1 kw  radiated  power  while  the  system 
described  employs  100  watts.  Combining  field  strength,  power,  antenna  gain, 
and  noise  yields  the  receiver  SNR. 

If  a receiving  antenna  with  high  gain  is  used  (a  certainty  if  military 
receiving  stations  are  used) , the  antenna  gain  can  be  added  to  the  S/N  ratios 
obtained. 

4.2.2  Ionospheric  Wave  HF  Telemetry  Link 

Ionospheric  waves  are  capable  of  propagation  over  much  longer  distances. 
However,  condi tions  vary  greatly  with  time  of  day,  year,  and  sunspot  cycle. 
Additionally,  many  modes  of  propagation  exist  in  the  high  frequency  band.  Con- 
sequently, performance  is  much  more  difficult  to  characterize  than  wdth  ground- 
waves. 
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(from  Picqucnard,  p.  220) 


Typical  Noise  Chart  for  HF 
Radio  Waves  on  Sumer  Day 
Between  0800  and  1200  Ilrs 
(from  Picquenard  p.  202) 
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Ionospheric  Modes 


The  primary  nodes  for  high  frequency  ionospheric  propagation  are  shown 
in  Figure  10,  with  the  main  ionospheric  layer  (i.e.,  F2~layer)  existing  at  a 
height  of  approximately  300  km. 

Ionospheric  Absorption 

This  is  one  of  the  major  factors  in  determining  signal  strength  over  a 
given  path.  It  is  characterized  by  the  absorption  coefficient  A,  which  is 
dependent  on  zenith  angle  of  the  sun,  season  of  the  year,. and  solar  activity. 
The  absorption  index  for  the  Gulf  Stream  area  of  interest  is  shown  in  Figure  11 
for  June  and  December  at  midnight  and  noon  for  extremes  of  the  sunspot  cycle. 
Figure  12,  from  Piequenard  (1974) , shows  how  the  signal  field  is  dependent  on 
the  absorption  index.  At  4.1  mHz,  the  signal  field  drops  rapidly  even  for 
small  values  of  A.  Therefore  ionospheric  propagation  over  distances  greater 
than  those  achievable  by  ground-v.ave  propagation  is  only  practical  at  night. 
Figure  13  shows  field  strength  as  a function  of  distance  frem  the  transmitter 
for  the  F2  skywave  mode  at  night. 

Availability  of  Skywave  Modes 

The  E and  F^  layers  are  not  useful  for  our  requirements.  The  equations 
for  E and  F^  layer  critical  frequencies  (i.e. , frequencies  above  which  a verti- 
cally-incident  wave  is  not  reflected)  are  given  as  follows : 

f E = 0.9  [(180  + 1.44  R)  cos  x)0,25 

f F.  = (4.3  + 0.01R)  cos0,20  X 
o 1 

where: 

x = zenith  angle  of  the  sun 

R = sunspot  number  (unimportant  for  this  writing) 

These  equations  show  that  the  critical  E and  F^  frequencies  approach  zero  when 
the  zenith  angle  readies  90°  (around  sunrise  or  sunset) . Therefore,  when 
absorption  decreases  enough  to  permit  propagation,  the  E and  F^  layers  have 
disappeared. 

The  F2  layer  behaves  well  for  this  application.  Propagation  should  be 
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Figure  10 


I IF  Skywave  propagation  Modes 
(from  Picquenard  p.  133) 


Month 

Tine 

x° 

Solar 

Zenith  Angle 

R12 

Sunspot 

Number 

A 

Absorption 

Index 

June 

0000 

- 

20 

0 

June 

0000 

- 

150 

0 

June 

1200 

15° 

20 

1.04 

June 

1200 

15° 

150 

1.50 

Dec 

0000 

- 

20 

0 

Dec 

0000 

- 

150 

0 

Dec 

1200 

60° 

20 

0.72 

Dec 

1200 

60° 

150 

1.13 

Figure  11  - Maximum  and  Minimum  Values  of  Absorption  Index 
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dB  above  1 


500  1000  1500  2000 

Distance  (km) 

Figure  13  - Signal  field  as  a Function  of  Distance 
for  1 kw  Radiated  Power 


possible  at  the  extemes  of  the  sunspot  cycle.  When  the  number  of  sunspots 
is  high,  the  critical  frequency  exceeds  4 mHz  throughout  the  evening,  as 
shown  in  Figure  14 . As  the  Itaximum  Useable  Frequency  (MUF)  will  always 
exceed  the  critical  frequency  (f  ) through  the  relation: 

MUF  = f sec  <}>  (4) 

c 

(0  = angle  of  incidence  at  the  reflecting  layer) 

therefore  the  mode  should  be  present. 

When  the  sunspot  number  is  lew,  the  critical  frequency  is  lower.  How- 
ever, it  is  still  adequate  for  propagation  at  4 mHz. 

In  order  to  explore  the  effect  of  oblique  incidence  on  MUF,  the  hypo- 
thetical case  shown  in  Figure  15 is  examined.  Ignoring  spherical  effects: 

" COS  * ~ [{D/2)2  + h2]*4  (5) 

For  a distance  of  600  km  with  an  Y^  layer  height  of  300  km  equation 
indicates  that: 

MUF  = 1.41  f 

Along  with  Figure  14  this  predicts  that  propagation  will  be  possible  through- 
out most  of  the  evening  (that  is,  when  the  critical  frequency  exceeds  3.0  mHz). 

When  the  buoy  is  within  about  200  to  300  km  of  the  receiving  station 
a different  situation  arises  as  indicated  in  Jenkins  et  al  (1967) . The  MUF 
for  ionospheric  propagation  is  not  raised  sufficiently  for  reliable  service 
and  hence  reliance  must  be  placed  on  the  ground-wave. 

Combined  Results 

Expected  signal-to-noise  ratios  throughout  the  range  of  Gulf  Stream 
interest  are  shown  in  Figure  16 . The  available  receiver  signal-to-noise 
ratios  are  computed  for  the  worst  time  of  year  at  the  best  time  of  day  (noon 
for  ground-waves,  evening  for  skywaves) . 
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Figure  14  - Critical  Frequency  (MHz)  for  F2-Layer 
Ionospheric  Reflection 
(from  Pioquonard  p.  124) 


Figure  15  - Definition  of  Incident  Angle  $ for 
One-Hop  F2 -Layer  Propagation 


HF  Receiver  Stations 


The  receiving  stations  indicated  in  Figure  1C  are  assumed  to  exist  at 
Canaveral, Florida; Cape  Hatteras,  Nantucket,  Mass. , and  Bermuda.  It  has  tenta- 
tively been  assumed  that  Navy  or  Coast  Guard  bases  could  be  errployed  on  a non- 
interference basis  for  data  acquisition.  In  such  cases  all  that  would  be  re- 
quired would  be  to  build  low  cost  automatic  data  acquisition  systems  that  turn 
on  every  12-hours.  The  basic  IF  signal  could  be  reinse'-ted  with  the  telemetered 
SSB  suppressed  carrier  signal  and  the  FSK  tone  put  on  a relatively  inexpensive 
tape  recorder  along  with  time.  The  tapes  could  then  be  converted  to  a com- 
puter-compatible format  at  the  laboratory  where  all  tapes  are  assembled. 

Because  a base  station  receiver  should  be  able  to  derive  useful  data 
from  as  low  as  a 10  db  signal,  or  lower , it  appears  from  Figure  16  that  the  full 
area  of  the  Gulf  Stream  shown  can  still  be  tracked  with  the  Florida,  Hatteras, 
and  Bermuda  stations.  If  only  the  Hatteras  and  Bermuda  stations  are  maintained 
the  HF  signal  may  be  quite  weak  in  the  region  of  the  Florida  Straits,  but  after- 
wards very  reasonable.  It  is  not  expected,  however,  that  the  type  of  higher 
frequency  trajectory  variability,  previously  found  in  a western  boundary  current 
such  as  the  Gulf  Stream,  would  exhibit  XwSclf  in  the  ccnfmod  region  of  the 
straits.  Therefore,  it  is  reasonable  to  initially  recommend  a two  base  station 
network  with  stations  at  Cape  -Hatteras  (possibly  at  Diamond  Light  shoals  tower 
^20  miles  offshore)  and  at  Bermuda  (possibly  Tudor  Hill) . 

The  cost  of  building  the  type  of  base  station  data  acquisition  system 
described  would  be  relatively  modest  (w$1000.-  replication  cost)  as  long  as  a 
receiver  and  a reasonable  receiving  antenna  are  available.  If  such  is  not  the 
case,  costs  per  station  would  increase  by  approximately  $4C00-$12,000.-  (com- 
ponent costs)  for  a receiver  and  antenna.  Appendix  3 describes  a few  pieces  of 
HF  equipment  that  enable  ball  park  system  cost  estimates. 

Buoy  Identification  Coding 

In  order  to  recognize  which  buoy  is  reporting  from  a batch  of  several 
buoys  deployed,  each  buoy  must  enploy  one  of  ten  different  300  Hz  bands  within 
the  3500  Hz  total  bandwidth  of  the  4.1625  to  4.1666  mHz  assigned  oceanographic 
band  (see  Livingston  et  al,  1977).  The  World  Administrative  Radio  Conference 
G'JARC)  has  assigned  this  bandwidth  which  allows  a 250  Hz  buffer  band  on  each 
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end  of  the  total  bond,  "here fore,  a maximum  of  ten  uniquely  identifiable 
buoys  working  this  band  can  be  in  the  water  at  one  time.  If  additional  buoys 
are  required,  the  6.2445  mHz  band  or  higher  must  be  used  and  depend  almost 
strictly  on  skywave  propagation. 

If  the  Teledyne  701  LORAN-C  receiver  with  memory  and  FSK  options 
is  used,  a standard  telemetry  feature  would  allow  for  8 identifications  bits 
such  that  each  HF  transmitter  could  work  at  the  same  frequency  if  desired. 

General  Comments 

It  should  be  reemphasized  that  the  IIF  system  suggested  is  somewhat  costly 
in  terms  of  development  funds  because  of  the  memory  required . In  order  to  in- 
crease overall  system  reliability  per  transmission  with  a single-frequency  trans- 
mission format,  the  candidate  IF  system  has  been  programmed  to  transmit  at  noon 
(using  groundwave)  ana  midnight  (using  skywave)  only.  With  such  a timeline, 
calculations  indicate  that  at  least  one  of  the  transmissions  will  be  reliable 
over  the  Gulf  Stream  range  of  interest.  An  alternative  approach  could  be  assumed 
in  which  periodic  LQPAN-C  TD's  were  immediately  radioed  to  si  ore  - requiring 
less  complicated  memory  and  timing  on  the  buoy.  This  approach  would  possibly 
require  the  transmittal  of  each  of  the  5 words  shown  in  Table  4 (i.e.,  52  bits) 
every  hour  at  a rate  of  100  bits/sec.  If  the  words  were  simply  assumed  to  be 
transmitted  once  per  hour  (for  example)  the  transmit  format  would  be  as  follows: 

HF  Hourly  Transmit  Format  (Mo  Memory) 


(1) 

Warm-up 

5.0 

sec 

(2) 

Carrier  Enable  + bit  sync 

7.5 

sec 

(3) 

52  bits/100  bits/sec 

.5 

sec 

13.0  sec 

If  done  hourly  at  a 200-watt  level  for  30  days  the  above  transmit  format 
would  lead  to  the  following  prohibitive  energy  consumption: 

JE^  (hourly  transmit)  = 1560  watt-hours 

Even  if  the  5-second  transmitter  warm-up  period  were  somehow  eliminated  the 
energy  consumed  would  still  be  approximately  960  watt -hours.  Therefore,  tlie 
data  format  described  in  Table  4 appears  to  be  the  most  energy,  and  by  all 
estimates,  the  most  oost-effective  way  to  go. 


Availability  of  Equipment 

Transmitter  costs  will  be  much  lower  for  a high  frequency  system  than 
for  a sattelite  system.  A large  number  of  mobile  transmitters  and  amplifiers 
have  been  developed  for  amateur  radio  use.  These  systems  are  reliable, 
require  a 12  vdc  power  supply,  and  cost  less  than  $400.  (see  Appendix  B) . 
Contact  has  been  made  with  TPL  Communications  (see  sheet  in  Appendix  B) , 
who  produce  suitable  power  amplifier  equipment  for  about  $250.  The  amplifier 
would  accept  a low  power  HF  signal  from  the  modulated  output  of  an  oscillator. 
The  total  parts  costs  for  an  HF  buoy  transmitter  system  (exclusive  of  timing, 
memory,  and  antenna  shown  also  in  Figure  7)  are  approximately  the  following: 


Transmitter  (i.e.,  oscillator  & modulator):  $400.- 
Power  Supply  $ 75.- 
Power  Amplifier  $250.- 


Tbtal  Estimated  Cost:  $725.- 


This  total  cost  can  be  compared  against  a $1500.  cost  for  a satellite  BTT,  but 
more  dollars  in  development  and  system  integration  costs  are  associated  with 
the  buoy  HF  transmitter  and  base  station  receiver.  Because  there  will  likely 
be  costs  associated  with  either  a daily  usage  ci large  for  the  TIEQS-N  satellite 
BTT  (i.e. , n,  $20/buoy  day)  or  the  purchase  of  a dedicated  TIROS-N  ground 
station  receiver,  the  costs  of  building  KF  transmitters  and  operating  two 
receiving  stations  are  not  overly  burdensome  by  comparison. 

Buoy  Antenna 

The  nost  practical  antenna  for  use  on  the  buoy  is  the  whip.  Such  an 
antenna  would  be  typically  a 2-3  meter  in  length  with  induction  loading  at  the 
base.  The  inductive  loading  is  necessary  to  make  the  8 ft  whip  resonate  at 
4.1  mHz.  For  a 4.1  mHz  transmission  frequency,  a vertical  antenna  would  have 
to  be  20-meter  high  if  no  loading  coil  was  used. 

In  this  application,  it  is  necessary  to  use  an  antenna  with  gain  that 
is  as  high  as  possible.  Equally  important,  the  antenna  must  not  lose  efficiency 
when  subjected  to  buoy  bobbing  and  rolling  motions. 
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For  an  inductively-loaded  whip,  high  efficiency  while  vertical  re- 
sults in  poor  efficiency  when  the  antenna  is  inclined  to  vertical.  Decreased 
efficiency  when  tilted  occurs  because  the  angle  change  with  respect  to  tiie 
conducting  plane  of  sea  water  results  in  a change  in  the  resonant  frequency 
of  the  antenna.  Unfortunately,  as  the  efficiency  of  the  antenna  is  increased, 
the  bandwidth  decreases.  Therefore,  the  change  in  resonant  frequency  wiiich 
occurs  when  the  antenna  is  tilted  degrades  the  performance  of  a high  effi- 
ciency antenna  more  titan  an  antenna  with  lower  efficiency. 

A theoretical  relationship  can  be  developed  by  wiiich  it  can  be  shown 
how  an  antenna  gain  varies  as  a function  of  Q for  various  buoy  tilt  angles. 

If  an  antenna,  with  Q = 50,  were  placed  on  a buoy  such  as  the  Scripps/McNally 
damped  spar  (see  McNally,  1976)  the  loss  of  antenna  gain  due  to  buoy  motion  is 
expected  to  be  less  than  3 db  in  a sea  state  5.  Such  would  be  true  because 
the  maximum  tilt  angle  is  expected  to  be  below  20  degrees.  Therefore,  by 
selecting  a 2.5-meter  whip  antenna,  tuned  to  a Q of  50,  the  total  effects  due 
to  buoy  motion  and  antenna  bending  are  expected  to  degrade  its  efficiency  by 
less  than  3 db.  Such  an  antenna  has  a gain  of  -15  db  when  upright.  Fbr  pur- 
poses of  analysis  a buoy  antenna  gain  of  -15  db  was  assumed. 

4.2.3  Skywave-Groundwave  Interaction 

Often,  signals  may -be  received  via  more  than  one  mode.  As  the  radio 
wave  trawls  different  distances  by  different  modes,  the  taro  signals  will  be 
received  at  slightly  different  times.  If  this  delay  approaches  the  time  dura- 
tion of  a single  bit,  destructive  interference  may  result.  As  a worst  case, 
interaction  between  a groundwave  signed,  and  the  same  signal  propagated  by  a 
double  hop  off  of  the  F^-layer  v/ill  be  considered.  If  the  signal  originates 
at  the  edge  of  the  groundwave  range  ( w 800  km) , assuming  an  F^-layer  height  of 
300  km,  the  time  delay  would  be  2.1  milliseconds.  If  data  is  transmitted  at 
100  bits  per  second,  or  1 bit  every  10  milliseconds,  the  skvwuve  would  arrive 
late  by  21%  of  a bit  cycle  period.  This  appears  to  represent  no  serious  problem. 
If  the  FSK  bit  rate  were  increased  the  problem  would  become  more  acute. 
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5.0  SUMMARY  AND  RECOMMENDED  SYSTEM  CONFIGURATION 


Table  5 presents  a summary  cost  estimate  of  five  possible  drifting 
buoy  system  combinations  which  should  meet  the  requirements  of  tracking 
the  higher  frequency  meanders  of  a western  boundary  current. 

Position  Determination 

The  study  has  shown  that  at  present  the  only  viable  means  of  acquir- 
ing buoy  position  data  at  the  desired  frequency  and  accuracy'  is  with  an 
automatic  DORAN-C  receiver.  It  has  been  tentatively  judged  that  the  Micro- 
logic ML-200  receiver,  although  possibly  not  of  the  sophistication  and  re- 
liability of  the  Northstar  6000  (from  Digital  Marine  Electronics,  Inc) , is 
adequate  for  the  job.  It  lowers  price  with  an  associated  penalty  of  a 
few  dollars  worth  of  extra  batteries  seems  like  a reasonable  choice.  This 
is  the  baseline  LORAN-C  receiver  chosen.  For  comparison.  Table  5 indicates 
the  use  of  the  Northstar  6000  (config.  3) . It  is  recosnnendcd  that  the  sys- 
tem schematically  shown  in  Figure  7 be  implemented  as  the  general  design. 
The  choice  of  a specific  LORAN-C  receiver  would  have  to  be  made  with  all 
of  the  most  recent  test  and  evaluation  data  at  Uut  time.  Tliis  is  true 
because  of  the  more  recent  construction  and  lack  of  test  data  on  models 
such  as  the  Northstar  6000  and  the  Intcrnav  LC-123 . If  the  choice  were 
to  be  made  with  the  data,  presented  herein  it  would  be  a difficult  decision 
between  the  Micrologic  ML-200,  because  of  its  low  price,  and  the  Teledyne 
701,  because  of  its  widespread  use  and  tests  by  the  Coast  Guard  (Cassis 
and  Adams,  1977)  plus  the  low  cost  options  which  include  a random  access 
memory  (RAM)  and  FSK  output,  which  have  already  been  implemented  for  the 
Department  of  Transportation,  Automatic  Vehicle  Monitoring  program  in 
Philadelphia  (Chambers  and  Stapleton,  1974).  It  is  felt  tlxat  the  choice 
at  a given  point  in  time  would  depend  on  a very  close  scrutiny  of  available 
test  data  combined,  possible,  with  a few  specific  tests  directly  pertinent 
to  performance  on  the  buoy  with  a TIROS-N  or  HF  transmitter  alongside. 

On-Board  Sensors 

The  other  main  cost  variations  in  Table  5 are  whether  or  not  a baro- 
metric pressure  sensor  is  to  be  used  on  the  buoy  and  how  the  data  are  tele- 
metered. Configuration  1 assumes  the  use  of  the  Bunker-Rano  sensor  being 
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Table  5 - Drifting  Buoy  Positioning  and  Data  Acquisition  System  Cost  Sunmary 


1 


employed  by  the  National  Center  for  Atmospheric  Research  (NCAR)  at  cost  of 
approximately  $500  each.  It  is  recctrmended  tint,  in  order  to  hold  down 
oosts,  the  LORAN-C/TIROS-N  system  eliminate  the  barometric  pressure  sensor 
and  use  configuration  2.  Because  of  their  relatively  low  cost  (^$100  each), 
it  is  assumed  that  both  a water  temperature  sensor  and  drogue  indicator  be 
included  on  the  buoy. 

Data  and  Position  Telemetry 

The  most  cost-effective,  reliable  manner  by  which  to  telemeter  the 
data  to  shore  appears  to  be  with  TIROS-N  satellite  transmitter.  If  a user 
cost  of  $20  per  buoy  per  day  is  imposed,  (at  6 buoys  installed  every  3 months 
for  a 30-day  lifetime  each  = 720  buoy-days  per  year)  an  additional  cost  of 
$14,400  per  year  is  expected.  This  represents  an  additional  $28,800  system 
cost  for  a 2-year  monitoring  program.  If,  on  the  other  hand,  am  automatic 
TIROS-N  ground  station  readout  system,  were  purchased  at  a cost  of  approxi- 
mately $20,000,  more  system  flexibility  would  be  achieved.  Such  a system 
would  require  programming  for  data  readout,  but  the  system  costs  could  po- 
tentially be  shared  by  other  experimenters  and  used  on  other  programs . 

This  approach  is  reconrnended . 

Configurations  4 and  5 assumed  data  telemetry  by  GOES  and  by  HF  re- 
spectively. The  GOES  system  is  deemed  too  costly  at  present.  The  LORAN-C/ 

HF  system  suggests  the  use  of  the  Teledyne  701  receiver  because  of  the  me- 
mory and  FSK  output  which  has  already  been  used  by  the  Department  of  Trans- 
portation. This  combination  can  be  purchased  for  approximately  $3,000.  The 
base  stations,  on  shore,  if  not  already  available  frem  the  Coast  Guard  or 
Navy,  could  be  set  up  by  erecting  a reasonable  antenna  (v$2, 000);  purchas- 
ing a good  HF  receiver  (<$ 5,000),  a Teledyne  FSK  single  sideband  demodula- 
tor (v$5,000),  and  a recording  system  (v$l,000).  A minimum  of  tars  such  sta- 
tions would  need  to  be  built  at  a total  estimated  store-side  component  cost 
of  $20,000.  Such  a system  is  also  competitive  on  buoy  system  costs  and  al- 
lows for  constant  buoy  position  determination  and  possible  retreival  if  de- 
sired. On  the  wtole,  such  a systan  would  be  more  cumbersome  to  operate  and 
less  reliable  than  the  LORAN-C/TIRGS-N  system.  A summary  of  the  overall 
system  costs  for  the  major  competing  system  choices  is  shown  in  Table  6. 

It  includes  not  only  the  purchase  of  capital  equipment  which  is  non-expendable 
but  also  rough  engineering  design  and  service  cost  estimates.  The  overall 
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sumnary  is  based  on  the  deployment  of  6 buoys  per  quarter  for  2 years. 

This  Table  does  not  include  overall  system  engineering  costs,  nor  does  it 
include  the  cost  of  the  fabrication  and  test  of  an  engineering  prototype, 
all  of  which  should  be  approximately  the  same  for  all  choices.  Table  6 
summarizes  only  a best  guess  at  the  incremental  costs  for  given  choices 
on  quantity  replication. 

Buoy  and  Drogue  Selection 

The  buoy  should  be  designed  such  that  the  antenna  systems  remain 
relatively  upright  (i.e.,  <±20°)  under  all  but  the  most  severe  conditions. 

This  is  not  to  eliminate  the  potential  use  of  a two-port  or  distributed 
buoyancy'  package,  both  of  which  might  result  in  better  drogue  performance  and 
survivability.  The  use  of  the  Scripps-McNa  1 lv  spar  may  be  a good  cost- 
effective  buoy  for  tlus  program.  This  specific  choice  should  be  nude  at  the 
time  of  the  design.  By  that  tiro  more  accurate  and  even  verified  math  models 
of  drogue  slippage  and  loads  may  be  available  to  help  in  the  choice. 

The  choice  of  a drogue  for  the  study  outlined  will  depend  very  much 
on  the  state  on  knowledge  of  drogue  performance  at  the  time.  If  no  more  is 
known  about  drogues  than  than  is  known  new,  it  is  tentatively  rcconmended 
that  bi-planar  crossed  vane  bo  employed  for  droguing  at  30  to  50-meter  depths . 

Future  Considerations 

The  system  schematically  shown  in  Figure  17 could  potentially  be  made 
more  flexible  through  the  use  of  modularity.  Figure  17  is  a general  block 
diagram  of  a modular  drifting  buoy  positioning  system  that  incorporates  many 
of  tiie  types  of  positioning  and  data  telemetry  systems  described  herein.  It 
can  be  made  to  function  anywhere  in  the  world  in  an  expendable  or  non-expendable 
mode  to  varying  degrees  of  accuracy  by'  interchanging  components.  For  example, 
Figure  17  encompasses  the  expendable  system  for  the  measurement  of  higher  fre- 
quency changes  in  western  boundary  currents.  By  changing  to  ■ n OMEGA  receiver 
and  an  HP  transmitter  an  experimenter  is  capable  of  monitori  j large  scale 
trajectory  data  anywhere  in  the  world  in  a non-expendable  mode  as  long  as  a 
radio  receiver  can  be  installed  within  approximately  1000-1500  km  of  the  ex- 
periment sight. 


Figure  17  - Modular  Concept  of  a Drifting  Buoy  Positioning  System 


6.0  ESTIMATED  BOUNDS  CN  DROGUED  BUOY  SLIPPAGE  ERRORS 


This  section  will  examine  some  of  the  current  basic  studies  of  the 
slippage  of  drogued  drifters  in  the  ocean.  In  addition,  current  research 
will  be  reviewed  which  permits  the  estimation  of  the  wind  and  water  drift 
forces  on  a buoy  - the  sum  of  which  gives  rise  to  the  error-inducing  slippage 
forces  on  a buoy.  The  wind  force  as  a percent  of  the  total  error  force  will 
vary  with  sea  state,  buoy  design,  wind  velocity,  surface  current  velocity, 
shear,  drogue  depth,  drogue  size,  and  drogue  effectiveness.  With  so  many 
parameters,  the  problem  is  very  corrplex.  Only  restrictive  assumptions  and 
bounds  can  be  put  on  the  problem  to  make  it  somewhat  manageable. 

It  wd.ll  be  assumed  for  analysis  purposes  that  there  is  no  current  at 
the  drogue  depth  and  that  the  drogue  is  shallow  enough  (V30-50  m.  depth)  and 
tethered  with  a small  enough  line  such  that  tether  line  horizontal  drag  forces 
can  be  neglected.  It  will  be  assumed  for  simplicity  that  the  wind  velocity 
can  be  readily  determined  by  a single  measurement  at  any  height  and  related  to 
that  at  the  buoy  by  the  Prandtl-von  Karman  universal  velocity  distribution. 

The  presence  of  the  air-water  interface  and  waves  are  assumed  to  produce  no 
disturbance  to  the  wind  field  around  the  buoy  such  that  the  buoy-wind  drag  co- 
efficients are  those  of  standard  cylinders  for  which  much  data  is  available  in 
Hoemer  (1965) . The  lateral  dimensions  of  the  buoy  are  further  assumed  to  be 
those  of  round  cylinders -of  diameter  much  less  than  a wavelength  (i.e.,  D«A) . 
the  NDBO  Nova  minibuoy  design  will  be  chosen  for  analysis.  The  assumed  base- 
line will  be  the  window  shade  supplied  by  NDBO  which  measures  approximately 
2.25  m.  x 9.83  m.  with  an  area  of  22  m2  and  a drag  coefficient  of  2.0. 


6.1  Buoy  Wind  Drag  Forces 

It  has  been  found  by  Vachon  (1975)  and  Saunders  (1976)  that  the  predom- 
inant error  force  in  most  drifting  buoy  applications  is  that  of  wind.  It  is 
reasonable  to  assume  that  the  rigid  above-water  portions  of  buoys  do  not  os- 
cillate in  a wind  field  in  the  manner  of  a mooring  line.  Iherefore  it  is  felt 
that  the  standard  round  cylinder,  open  flow  drag  coefficient  of  approximately 

1.1  to  1.2  applies  fairly  well  at  subcritical  Reynolds  numbers  (cylindrical 
mooring  lines  may  strum  leading  to  much  higher  effective  drag  coefficients) . 
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If  the  ratio  of  buoy  height-to -diameter  is  of  the  order  of  5 or  less  it 
appears  that  the  cylinder  drag  coefficient  may  be  as  lew  as  .73  at  sub- 
critical  Reynolds  numbers  (Hocmer,  1965,  p.  4-3).  In  spite  of  these  data, 
values  of  CD  = 1.0  and  1.1  have  been  used  by  Vachon  (1975)  and  CD  = 1.1  by 
Saunders  (1976)  for  portions  of  the  Nova  buoy  subject  bo  wind  forces.  These 
assunptions  will  be  used  in  developing  a maximum  estimate  of  slippage.  The 
analysis  by  Vachon  (1975)  found  that  either  the  values  used  for  wind  drag 
were  too  lav  or  the  value  of  the  drogue  drag  coefficient  was  too  high.  The 
latter  was  generally  concluded. 

CASE  I - Wind  Force  Only  - No  Shear 

If  an  experimenter  were  to  estimate  drifter  trajectory  deviations 
caused  only  by  wind  influences  on  the  non-wetted  portion  of  a buoy,  an  inter- 
mediate bound  can  be  placed  on  the  errors.  This  case  might  well  be  the  only 
type  of  correction  that  could  be  made  based  on  avai  lable  wind  field  data  and  a 
lack  of  knowledge  of  the  surface  current,  wave  influence,  and  drogue  condition 
or  behavior.  For  estimation  purposes  the  Nova  minibuoy  outlined  in  Figure  18 
will  be  assumed  with  the  estimated  cross-sectional  areas  and  drag  areas  shown. 
Figure  18  also  indicates  the  Scripps-McNally  spar  with  estimated  cross- 
secitonal  areas  and  drag  areas  given  for  comparison.  These  buoys  will  exper- 
ience tlie  estimated  wind-  forces  given  in  Table  7. 

With  no  shear  assumed  between  the  surface  and  drogue  depth  it  is  possible 
to  sum  the  drag  areas  of  the  wetted  portion  of  the  buoy  and  the  drogue  (i.e., 
(CpA)  = 44.74  m2).  Both  wetted  elements  of  the  system  then  act  as  a drogue, 

TABLE  7 

Estimated  'Wind  Forces  on  Drifting  Buoys 
Wind  Vel.  Nova  Forces  Scripps/McNally  Spar 


knots 

ny/s 

(Newtons) 

lbs 

(Newtons) 

lbs 

10 

5.1 

7.4 

1.7 

5.5 

1.2 

20 

10.3 

30.2 

6.8 

22.3 

5.0 

30 

15.4 

67.6 

12.9 

49.8 

11.2 
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Figure  18  - Estimated  Cross-Sectional  and  Drag  Areas 

of  Drifting  Buoys 
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restricting  the  error-inducing  wind  force.  Ihe  resultant  estimated  system 
slip  is  shown  in  Table  8 for  the  Nova  minibuoy  only.  The  slip  of  the  Scripps 
buoy  can  be  calculated  by  using  values  shown  in  Table  7 . 

TABLE  8 

Estimated  Nova  Buoy  Slip  with  No  Shear  Present,  (CL)  = 2.0 

’ ~ drogue 

Wind  Vel.  Slip  Max.  Trajectory  Error/hr 


(knots) (M/s) (n\/s) (m) 


10 

5.1 

.018 

±65 

20 

10.3 

.036 

±130 

30 

15.4 

.054 

±195 

If  the  drogue  is  working  ineffectively  because  it  is  either  not  weathervening 
proper lv,  in  the  case  of  a window  shade  drogue/  or  is  hanging  downward  or  torn, 
in  the  case  of  a parachute,  then  its  drag  coefficient  averaged  over  time  may 
drop  to  approximately  0.5.  This  condition  would  lead  to  the  slippage  values 
shown  in  Table  9.  It  should  be  noted  that  the  drag  coefficient  assumed  has  a 
basis  in  measured  data.  In  Vachon  (1973) , the  drag  coefficient  of  a porous 
fishing  net  was  evaluated  to  be  approximately  . 2 to  .3  by  towing  it  parallel 
to  the  flow  direction.  It  is  assumed  that  the  window  shade  would  weathervane 
properly  some  of  the  time  - raising  the  average  as  indicated.  In  addition,  a 
parachute  can  still  function,  but  at  a lower  level  of  effectiveness,  even  though 
hanging  downward  or  tom. 

TABLE  9 

Estimated  Nova  Buoy  Slip  with  No  Shear,  (CL)  = .5 

drogue 

Wind  Vel.  Slip  Max.  Trajectory  Error/hr 


(knots) 

(mV's) 

(m/s) 

(m) 

10 

5.1 

.036 

±129 

20 

10.3 

.072 

±261 

30 

15.4 

.108 

±390 
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It  can  be  seen  by  comparing  Tables  7 and  8 that  the  drogue  could 
function  ineffectively  in  moderate  winds  and  the  drifter  trajectory  errors 
will  only  begin  to  exceed  the  desired  accuracy  of  the  LORAN-C  positioning 
system  with  30-knot  winds.  Even  with  an  effective  drogue  the  trajectory 
errors  can  became  large  with  over  30-knot  winds. 

6.2  Drag  on  Wetted  Portions  of  a Buoy 

The  forces  on  the  wetted  portion  of  the  buoys  described,  although 
generally  estimated  to  be  of  smaller  order  than  the  wind  forces,  are  felt  to 
be  much  higher  than  simple  calculations  indicate  (Vachon,  1975  and  .Saunders, 
1976). 

CASE  II  - Wind  Driven  Surface  layer 

Wu  (1975)  has  considerable  data  to  suprxsrt  a claim  that  for  the  case  of 
no  major  current  regimes,  such  as  the  Gulf  Stream  , the  steady  wind- induced 
surface  current  will  equal  3.5  percent  of  the  wind  velocity  independent  of 
fetch,  at  fetch.es  greater  than  approximately  100  km.  For  shorter  fetches  the 
ratio  will  increase  to  approximately  4£  at  wind  velocities  of  20  m/s.  Wave- 
induced  Stokes  drift  current  was  estimated  to  constitute  approximately  one- 
seventh  of  this  (i.e.,  . 5 p:).  For  the  analysis  herein,  it  will  be  assumed  that 
this  current  decays  as  e ( k = ^/\)  such  that  the  current  forces  on  the 
buoy  can  be  estimated.  For  the  purposes  of  analysis  the  Nova  minibuoy  will  be 
divided  into  three  wetted  parts  as  described  in  Table  10.  As  the  wind  velocity 
is  varied  the  wavelength  of  the  seas  is  assumed  to  vary  in  accordance  with  the 


TABLE  10 


Nova  Minibuoy  Estimated  Wetted  Drag  Characteristics 


Buoy 

Section 

Cross-Sectional 
Area,  A (m2) 

Z.  = Depth 

1 ' (m) 

°D 

cda 

(m2) 

K.  = ^(CdA)  j 
(Kg  • m *)  I 

(1)  Cone 

.3 

• 

CJ 

CO 

1 

.24 

123.2 

(2)  Narrow  Cyl. 

.23 

1.37 

KH 

.25 

128.4 

(3)  Battery 

.14 

2.44 

CSS 

.13 

66.7 

Totals 



.62 

318.3 

k i 


Pierson-toskowitz  spectrum.  For  the  three  wind  velocities  previously  exam- 
ined, Table  11  presents  a summary  calculation  of  the  wind-driven  current  forces 


TABLE  11 

Current  Forces  on  Nova  Buoy 


Buoy 

Section 

Wind  Velocity 

(knots)  OVs) 

Surf. 

Curr. , V 
(Vs) 

s 

> 

_ _ -kz. 

= V e l 
s 

(Vs) 

K.V2  (Z) 
l 

(Newtons) 

(1) 

10 

5.1 

.179 

35 

.17 

3.4 

(2) 

10 

5.1 

.179 

35 

.14 

2.5 

(3) 

10 

5.1 

.179 

35 

.12 

2.5 

(1) 

20 

10.3 

.361 

100 

.35 

15.1 

(2) 

20 

10.3 

.361 

100 

.33 

14.0 

(3) 

20 

10.3 

.361 

100 

.31 

6.4 

(1) 

30 

15.4 

.539 

300 

.53 

34.6 

(2) 

30 

15.4 

.539 

300 

.52 

34.7 

(3) 

30 

15.4 

.539 

300 

.51 

17.5 

Sumration  of  Current 

Drag  Forces 

(N) 

10  - 

knot  wind:  = 

6.8 

20  - 

knot  wind:  = 

35.5 

30  - 

knot  wind:  = 

86.8 

on  the  buoy.  The  sum  of  these  forces  over  the  three  buoy  sections  is  listed 
at  the  bottom  of  Table  11.  In  the  most  extreme  case,  these  forces  are  assumed 
to  always  act  col  incar  and  sum  with  the  wind  forces  given  in  Table  7 . Table  12 
presents  the  new  buoy  force  sum  as  well  as  the  expected  slip  and  trajectory 
deviation  for  a drogue  with  drag  coefficients  of  2.0  and  .5  . It  is  apparent 
in  this  analysis  that,  if  no  trajectory  corrections  are  made  based  on  wind  data, 
the  hourly  trajectory  error  will  be  of  the  same  order  as  the  LORAN-C  position 
error  in  the  presence  of  up  to  20  knot  winds.  Above  this  value,  wind  field  data 
must  be  acquired  and  corrections  made  in  order  to  hold  dewn  the  errors.  If 
the  drogue  is  assumed  to  be  working  properly  {i.e.,  ~ 2.0),  but  in  reality 
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TABLE  12 


Maximum  Buoy  Forces  and  Drogue  Slippage  for  Sum 
of  Wind  and  Wind-Induced  Current  Forces 


Wind  Velocity 
(knots)  (m/s) 

Force  Sim 
(Newtons ) 

Drogue  Slio 
(Cd=2.0) (m/s) 

Itajectory 
Error/hr  (m) 

Drogue  Slip 
(Cq=*5)  (Try's) 

Trajectory 
Error/hr  (m) 

10 

5.1 

14.2 

.024 

±90 

.05 

-179 

20 

10.3 

65.7 

.053 

±193 

.107 

±385 

30 

15.4 

154.4 

.082 

±295 

.164 

±590 

it  is  not  (i.e.,  CQ  ~ .5),  the  error  in  the  trajectory  correction  itself  will 
vary  from  89  m.  (at  10  knots)  to  295  m.  (at  30  knots).  For  this  analysis,  the 
error  in  drogue  drag  coefficient  gives  rise  to  trajectory  correction  errors 
equal  to  the  first  order  trajectory  correction  itself  if  the  drogue  was  working 
properly. 

The  v.orks  of  Keulegan  and  Carpenter  (1958)  and  Sarpkaya  (1976)  were 
employed  in  order  to  analyze  how  the  estimates  of  wetted  drag  coefficients  might 
vary  in  wave-induced  oscillatory  tier/;  compared  to  those  values  in  Table  10  derived 
under  steady  flaw  conditions.  This  analysis  lead  to  approximately  a 10T  in- 
crease in  CpA  because  the  C^  in  the  cylinder  area  was  raised  to  1.35  while  the 
others  were  unchanged.  The  additional  trajectory  correction  required  is  mini- 
mal. If,  however,  one  is  beginning  to  estimate  the  various  C^'  s for  wetted 
portions  of  a buoy  these  analyses  should  be  employed  rather  than  values  derived 
under  steady  flow  conditions. 

The  work  of  Saunders  (.1977)  could  also  be  employed  in  the  case  of  oscilla- 
tory flow  in  order  to  get  an  estimate  of  wetted  forces.  His  werk  applies  to 
forces  in  the  presence  of  non-zero  currents  for  which  |awj»  [U  | (a  = wave 
amplitude^/  ■=  frequency,  ju  | = average  surf,  current).  His  main  conclusion  is 
that  for  many  applications  the  average  drag  on  a wetted  body  in  oscillatory  flew 
my  be  written  with  sufficient  accuracy  by  the  expression: 

F = (hp  C^A)  aw  U (6) 

He  is  mainly  addressing  the  problem  of  averaging  a velocity  relative  to  the  buoy 
which  is  composed  of  a steady  and  an  oscillatory  flew  component  as  follows: 
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(7) 


V jjv  jJ  = [U  + a w cos  (kx  - wt)]|u  + aw  cos  (kx  - wt)  | 

This  relation  can  be  simply  averaged  by  integration  as  long  as  |u|>|aw|.  When 
such  is  not  the  case  it  must  be  done  numerically  or  use  equation  (3)  as  an 
approximation.  The  value  of  employed  in  (5)  should  come  from  the  work  of 
Keulegan  and  Carpenter  (1958)  and  Sarpkaya  (1976).  The  data  presented  in 
Tables  11  and  12  were  not  redone  using  equation  ( 6)  because  the  general  results 
of  this  cliapter  would  change  little. 

6.3  Comparison  with  r-ieasured  Slippage  Data 

Two  instrumented  slippage  experiments  have  been  conducted  on  drifting 
buoys  in  recent  years.  Both  Vachon  (1975)  and  Saunders  (1976)  have  indicated 
an  inability  to  get  closure  on  a force  balance  on  a drifting  buoy  unless  either 
the  drogue  is  of  a drastically  lower  drag  coefficient  tlian  derived  under  steady 
flow  tests  or  wave  forces  are  greater.  In  Vachon  (1975)  there  are  indications 
that  a window  shade  drogue  may  not  weathervane  properly  v.iiich  supports  the 
former  contention.  Furthermore , James  McCullough  at  Foods  Hole  (personal  con- 
tact) has  demonstrated  a strong  correlation  between  measured  slip  and  wind 
velocity.  In  both  experiments  insufficient  ground  truth  data  were  available  by 
which  to  ascertain  true  current  at  the  drogue  depth  or  true  shear  forces  on  the 
buoys.  As  a result  there  is  at  present  no  firm  slippage  data  which  is  fully 
correlatable  to  the  forcing  on  the  buoy. 

6.4  Drogue  Performance 

Because  of  present  questions  regarding  the  performance  of  both  parachutes 
and  window  shade  drogues,  it  seems  reasonable  to  mount  a bi -planar  crossed  vane 
drogue  to  a buoy  for  the  study  of  western  boundary  currents.  Tliis  drogue,  with 
a scale  model  drag  coefficient  of  C^  = 1.13  (Vachon,  1973)  should  behave  in  a 
more  predictable  fashion  allowing  adequate  trajectory  corrections  from  wind  field 
data.  As  more  positive  performance  data  are  derived  on  the  other  drogues  - they 
should  be  used  where  feasible. 
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APPENDIX  A 


TIROS-N  Ebrd  structure  and  Data  Rato  Description 


For  purposes  of  estimation  it  will  be  assumed  that  every  12  hours  the 
satellite  transmitter  dumps  the  previous  12-tours  of  LORAN-C  data.  At  these 
times  it  will  be  assumed  that  it  transmits  a 256-bit  ward  stream  for  the 
parameters  shown  in  Table  A-l.  In  order  to  code  the  LORAN-C  time-difference 


Table  A-l 

TIROS-N  Satellite  Transmitter  Data  Word  Structure 


Parameter 

Bits/Word 

Max.  Words/XMSN 

Total  Bits 

LORAN-C,  TD. 

12 

9 (i.e,  1 pt./80') 

108 

LORAN-C,  TD2 

12 

9 

108 

Water  Temperature 

12 

1 

12 

Barometric  Pressure 

12 

1 

12 

Drogue  Indicator 

4 

1 

4 

Total  useful  bits/12  hrs: 

224 

words  (TD's)  into  a 12-bit  format,  tire  first  two  digits  of  the  normal  6-digit 
array  (i.e.,  5-digits  plus  single  decimal  digit)  must  be  dropped.  In  other 
words  if  the  normal  TD  were  to  be  69,046.2  (microseconds),  the  first  two  digits 
could  be  easily  dropped  while  still  retaining  position  knowledge.  This  is 
possible  as  long  as  it  is  known  which  station  pair  is  being  tracked  and  there 
are  no  long  periods  during  wliich  the  buoy  position  is  not  approximately  known. 

The  244-bits  of  data  being  transmitted  every  12-tours  to  the  TIROS-N 
satellite  represents  a pair  of  buoy  TD’s  every  80  minutes.  The  standard 
TIROS-N  format  is  expected  to  be  32,  8-bit  words  (i.e.,  256  bits)  with  every 
transmission  The  extra  12  bits  capable  of  being  added  to  the  244-bit  format 
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in  Table  A-l  could  represent  a designation  of  LORAN  Station  pairs  or  another 
data  ward  such  as  battery  voltage  or  a LORAN-C  signal  irenitoring  function. 
Besides  the  256-bits  per  transmission,  the  buoy  satellite  transmitter  would 
send  up  additional  bits  for  synchronization  and  platform  identification.  As 
with  the  NIMBUS-6  satellite,  the  TIROS-N  buoy  transmitters  can  be  multiplexed 
such  that  the  first  256-bits  could  represent  one  set  of  data  and  approximately 
one  minute  later  the  next  256  bits  represents  the  remainder  of  the  data.  With 
such  a feature,  a full  24-hours  of  data  can  be  theoretically  stored  in  the 
suggested  1024  bit  random  access  memory  (RAM)  and  dumped  every  12-tours  - 
allowing  a double  redundancy  in  daca  transmission . 
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APPENDIX  B 


IF  Transmission  Equipment 

There  is  a limited  amount  of  IF  transmission  equipment  available  for 
the  purposes  described  in  this  writing.  The  following  are  examples  of 
ham-type  equipment  that  might  be  employed  for  transmitting  data  from  the  buoy 
and  conmarical-type  equipment  that  could  be  purchased  for  the  receiving  station. 

Buoy  IF  Transmitter  Equipment 

(1)  TEN-TEC  Triton  IV  transceiver 

- 200  watts  input 

- freq.  3.5  to  29.7  MHz 

- SSB  with  12  VDC  drive  9 $699. 

(2)  TPL  Solid  State  Linear  amplifier 

- 1 watt- in,  70  watts-out. 

- SSB  @ v $170.- 

(3)  Heathkit  Model  IIW-104  SSB 

5 band  transceiver  (3.5  - 29  14Hz)  - $490.- 


Base  Station  Receiving  Equipment 


(1) 


IF  Receivers  - 

National  Radio  Corp. 
Collins  Radio 
Others 


s,  $5,000. 

s,  $8,000  - $11,000. 
% $8,000  - $12,000. 


(2)  Demodulators: 

Teledyne : 


$5,000. 
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